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ABSTRACT 


A  novel  high-speed  jet  noise  suppression  technique  using  high-pressure  gas 
microjet  injection  at  the  nozzle  exit  was  developed  with  promising  results  using  the 
laboratory  scale  jet.  The  main  jet  parameters,  the  nozzle  pressure  ratio  and  the 
temperature  ratio  are  chosen  to  correspond  with  realistic  engine  operating  conditions. 
Keeping  in  mind  the  applicability  of  the  technique  to  full-scale  engines,  the  microjet  mass 
flow  was  kept  at  less  than  2%  of  the  primary  jet  mass  flow.  The  A-weighted  spectrum  with 
appropriate  scaling  to  reflect  the  full-scale  nozzle  show  a  6  dBA  reduction  in  the  peak 
noise  radiation  direction.  An  attempt  is  made  during  this  investigation  to  develop  Active 
Noise  Control  methods  for  cancellation  of  broadband  jet  noise.  A  demonstration  of 
successful  control  architecture  has  been  accomplished  using  the  colored  random  noise 
cancellation.  An  actuator  based  on  pulsed  microjet  appears  to  be  most  promising  for  the 
current  application.  A  specially  designed  micro  valve  operating  at  high  pressures  and 
high  frequencies  was  developed  and  tested  successfully.  The  valve  was  able  to  provide  a 
pulsed jet  operating  at  nozzle  exit  mean  pressure  of  about  45  psia  at  2000  Hz. 
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1.  INTRODUCTION 

This  report  presents  the  results  of  active  control  techniques  used  to  suppress  the 
far-field  noise  emanated  from  jet  exhaust  of  fighter  aircraft,  such  as  the  F- 18  (figure  1),  at 
take-off  and  landing.  This  technology  was  developed  with  a  view  that  we  do  not  have 
easy  access  to  the  primaiy  sources  of  sound-  the  unsteady  pressure  fluctuations  in  the  jet. 
Instead  an  attempt  is  made  to  control  the  acoustic  fields  generated  by  the  primaiy  sources 
by  the  operation  of  secondary  sources.  The  control  will  be  over  the  near  field  of  the 
primaiy  sources  such  that  the  far  field  acoustic  power  is  reduced  or  that  the  directivity  of 
the  field  is  substantially  altered.  Destructive  interference  has  long  been  appreciated  as  the 
source  of  great  character  in  acoustic  fields.  “Knowing  that  two  different  source 
distributions  can  generate  the  same  wave  field  and  one  source  distribution  is  under  our 
control,  then  a  simple  change  in  sign  makes  the  primaiy  noise  field  subject  to  extinction 
by  the  presence  of  the  secondary.  Furthermore,  even  though  the  source  of  anti-noise  can 
be  of  completely  different  construction  to  that  of  the  primary  field,  the  silence  is,  in 
principle,  achievable  eveiywhere  outside  the  source  distribution  *. 


Figure  1.  F-18  aircraft  in  flight  with  visible  nozzle  exhausts. 

Active  noise  control  methods  have  been  successfully  applied  to  attenuate  the 
cabin  noise  in  propeller-driven  aircraft.  It  is  to  be  noted  that  the  cabin  noise  is  highly 
tonal  and  their  levels  decay  rapidly  with  increasing  mode  number.  As  a  result,  a  small 
number  of  sensors  and  actuators  are  required  to  achieve  significant  levels  (~10dB)  of 


*  J.  E.  Fffowcs  Williams,  “Anti-sound”,  Proc.  R.  Soc.  Lond.  A  395, 1984,  pp.  63-88. 
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noise  reduction  in  the  cabin.  More  recently,  the  control  of  combustion  instabilities  has 
received  considerable  attention  because  of  large  amplitude  flow  oscillation  in  the 
combustion  chamber  can  often  have  very  deleterious  effects.  In  these  cases  an  unstable 
feed  back  cycle  is  present,  with  part  of  the  cycle  mediated  by  pressure  oscillations  in  the 
combustion  chamber.  The  aim  of  the  active  control  here  is  to  change  the  effective 
behavior  of  the  cycle  to  achieve  stability. 

It  is  important  to  stress  that  the  active  control  is  most  feasible  and  appropriate  at 
low  frequencies,  where  the  near  field  is  extensive.  For  example,  a  typical  far  field 
microphone  signals  of  a  full-scale  engine,  operating  at  conditions  similar  to  that  of  an  F- 
18  engine  depicts  low  frequency  components  with  peak  Strouhal  number  (fD/U;  f: 
frequency,  D:  Nozzle  exit  diameter  and  U:  Nozzle  exit  velocity)  being  around  0.2.  The 
dominant  frequencies  are  within  the  range  of  100  ~  1000  Hz.  At  these  relatively  low 
frequencies,  conventional  passive  attenuation  schemes  involving  dissipation  of  sound 
energy  by  porous  or  microstructured  materials  become  impractical  and  active 
cancellation  techniques  offer  for  them  the  best  prospect  of  control.  However,  the  high¬ 
compactness  ratios  (=  k  l,  where  k  =  f/c,f  is  the  frequency  and  c  is  the  speed  of  sound 
and  l  is  a  typical  source  size)  of  the  sources  involved  require  carefully  designed  actuators 
and  their  implementation. 

In  this  work  we  attempted  to  inhibit  the  radiation  from  the  source  region  by  the 
use  of  microjet  injection  using  high-pressure  air  and  water.  During  this  period  attempts 
are  also  made  with  active  noise  control  scheme  using  high  frequency  actuators.  One  of 
the  technical  challenges  faced  during  this  investigation  was  the  development  of  a  high 
frequency  actuator  with  sufficient  power  to  affect  the  initial  jet  and  resulting  acoustic 
characteristics. 

The  goal  of  this  effort  is  to  demonstrate  10  dB  broadband  noise  reduction  at  the 
maximum  noise  radiation  angle  of  a  high  speed  supersonic  jet  operating  at  conditions 
similar  to  that  of  F-18E/F.  Next  section  describes  our  technical  approach  and  the  results. 
Next  section  describes  the  control  system  architecture,  including  the  several  actuator  and 
sensing  systems  and  the  template  for  the  adaptive  control  scheme. 
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2.  HIGH  SPEED  JET  NOISE  SUPPRESSION  WITH 
MICROJET  INJECTION 


The  far  field  noise  of  a  supersonic  jet  is  comprised  of  four  major  noise 
components1.  The  first  is  a  high  frequency,  short  wavelength  field  that  is  coherent  in 
phase,  commonly  referred  to  as  Mach  waves.  They  have  plane  phase  fronts  and  are 
confined  to  a  definite  wedge  sector  and  emanate  from  the  region  within  the  first  few 
diameters  downstream  of  the  nozzle  exit,  as  can  be  seen  in  Figure  2.  These  are  generated 
by  small-scale  disturbances  (or  eddies)  that  are  being  convected  at  supersonic  speeds  so 
that  they  emit  Mach  waves  in  the  direction  defined  by  a  disturbance  convection  velocity 
and  the  atmospheric  speed  of  sound2,3,4.  Surrounding  the  jet  with  a  gas  stream  that  has  a 
higher  speed  of  sound  eliminates  these  waves  as  demonstrated  by  Oertel  and  Patz5  and 
more  recently  by  Papamoschou6.  The  second  field  is  highly  directional,  peaking  at 
smaller  angles  to  the  jet  axis  (or  larger  angles  to  the  inlet  axis).  This  noise  field  is 
generated  from  large-scale  instabilities  reaching  peak  amplitudes  in  the  region  somewhat 
upstream  of  the  end  of  the  potential  core.  These  sources  of  noise  are  associated  with  the 
unsteady  flow  on  a  scale  that  is  comparable  with  the  local  shear  layer  width.  The  spectral 
intensity  of  this  sound  field  constitutes  two  distinct  peaks7.  One  is  associated  with  the 
highly  directional  Mach  waves  characterized  by  high  positive  pressure  peaks  in  the  far- 
field  microphone  signal8.  These  Mach  waves  are  of  significant  strength  as  compared  to 
those  that  originate  very  close  to  the  jet  exit  as  discussed  above.  This  intense  radiation  is 
observed  to  emanate  from  a  region  between  5~10  nozzle  diameters  and  is  associated  with 
supersonically  traveling  large-scale  coherent  regions  of  vorticity9  (Figure  1).  It  is  found 
that  the  far-field  intensity  contribution  of  this  source  is  about  30%  of  the  measured  total 
intensity7,10.  The  sources  for  the  second  peak  appear  to  be  located  much  further 
downstream  (10-20  nozzle  diameters)  and  are  associated  with  the  unsteady  flow 
generated  by  the  large  structures  quite  similar  to  that  in  subsonic  jets. 

The  third  noise  field  is  at  all  angles  to  the  jet  axis  and  at  higher  frequencies.  This 
sound  is  generated  in  precisely  the  same  manner  as  in  subsonic  flow  by  the  conventional . 
chaotic  turbulence.  Recent  analysis  of  experimental  data  by  Tam11  indeed  shows  vividly 
the  contributions  of  the  distinct  components  to  the  total  far  field  spectrum— that  associated 
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with  large-scale  motions  (inclusive  of  Mach  wave  radiation  if  it  exists)  and  a  component 
to  small-scale  turbulence. 

The  fourth  noise  field  is  commonly  referred  to  as  shock-associated  noise  and  it 
occurs  in  non-ideally  expanded  jets.  The  far-field  noise  spectrum  associated  with  this 
noise  source  typically  consists  of  discrete  peaks,  representing  the  screech  tones  and  a 
broad  peak  associated  with  the  shock-associated  broadband  noise12. 


Figure  2.  Schlieren  picture  of  a  Mach  2  round  jet.  To=  1250  K.  Uj  =  1050  m/s 
This  section  presents  a  possible  approach  for  the  suppression  of  the  dominant 
large  scale  mixing  noise  sources  in  a  supersonic  jet.  High-pressure  microjets  are  injected 
into  the  primary  jet  at  the  nozzle  exit  to  manipulate  the  dominant  source  region,  which 
extends  typically  from  5  to  20  diameters  from  the  nozzle  exit.  Recent  results  of 
experiments  on  an  Mj  =  0.9  round  jet  suggests  that  the  interaction  of  the  microjets  with 
the  jet  shear  layer  reduces  the  turbulence  levels  in  the  noise  producing  region  of  the  jet13. 
It  appears  that  the  microjets  influence  the  mean  velocity  profiles  such  that  the  peak 
normalized  vorticity  in  the  shear  layer  is  significantly  reduced  and  thus  inducing  an 
overall  stabilizing  effect.  Therefore,  it  is  suggested  that  that  an  alteration  in  the 
instability  characteristics  of  the  initial  shear  layer  can  influence  the  whole  jet  exhaust 
including  its  noise  field. 
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The  present  experiments  were  conducted  using  a  convergent  axisymmetric  nozzle 
operating  at  nozzle  pressure  ratios  of  3,  which  resulted  in  jet  Mach  numbers  of  1.38.  The 
nozzle  temperature  ratio  (stagnation  temperature/ambient  temperature)  was  kept 
nominally  at  3.  The  microjet  mass  flow  rate  was  kept  at  less  than  2%  of  the  main  jet 
mass  flow  rate. 

2. 1  Apparatus,  Instrumentation  and  Procedures 

Experiments  were  conducted  in  the  newly  built  High  Temperature  Supersonic  Jet 
Facility  at  the  Fluid  Mechanics  Research  Laboratory  of  the  Florida  State  University.  A 
schematic  of  the  facility  can  be  seen  in  Figure  3.  A  high-displacement  reciprocating  air 
compressor,  which  is  capable  of  supplying  air  at  a  maximum  storage  pressure  of  160 
bars,  drives  the  facility.  Large  storage  tanks  provide  a  total  capacity  of  10  m3.  After 
leaving  the  storage  tanks,  the  air  is  subjected  to  a  two-stage  pressure  control  system 
before  it  is  heated  by  passing  through  a  Sudden  Expansion  (SUE)  Burner  that  uses 
ethylene  as  the  fuel.  The  burner  was  originally  designed  to  supply  the  high  Mach  number 


Figure  3.  A  Schematic  of  the  High  Temperature  Supersonic  Jet  Facility. 
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flow  for  the  testing  of  ramjet  engines.  The  burner  has  a  simple  and  reliable  construction 
and  it  also  has  wide  flame  stability  limits  due  to  the  fact  that  the  fuel  injectors  are  located 
at  the  point  of  expansion.  The  burner  has  an  inlet  diameter  of  76.2  mm  and  an  exit 
diameter  of  152.4  mm  with  an  overall  length  of  0.813  m  and  it  is  cooled  by  a  water  jacket 
that  surrounds  it. 

2.  la  Flow  Management 

To  minimize  the  effect  of  any  combustion  instabilities  on  the  flow,  a  1.525  m  long 
water-cooled  pipe  having  a  diameter  of  152.4  mm  is  used.  Following  the  water-cooled 
extension  are  three  sections:  a  measurement  section,  an  extension  section  and  a  reduction 
section.  Each  of  them  is  made  of  a  nickel-based  alloy,  Nistele  230,  which  can  withstand 
temperatures  in  excess  of  1400  K  without  the  need  for  cooling.  The  conditions  of  the 
flow  before  reaching  the  nozzle  are  carefully  monitored  and  managed  within  these 
sections. 

Immediately  following  the  water-cooled  extension  is  the  measurement  section. 
This  section  is  0.61  m  long  and  has  an  inside  diameter  of  152.4  mm.  The  first  50  mm  of 
this  section  contains  a  ceramic  honeycomb  followed  by  two  titanium  screens  to  suppress 
any  large-scale  disturbances  in  the  flow.  Midway  into  the  section  are  four  equally  spaced 
ports  along  the  circumference  that  are  2.54  mm  in  diameter.  Two  of  the  opposing  ports 
are  used  to  measure  the  static  pressure  of  the  airflow.  It  turns  out  that  at  this  location  the 
static  and  stagnation  pressures  are  nearly  equal.  The  other  two  opposing  ports  have 
Omega  C-type  thermocouples  inserted  through  them  so  that  they  can  measure  the 
stagnation  temperature  of  the  airflow.  The  extension  section  follows  the  measurement 
section.  This  section  is  also  0.61  m  long  and  has  an  inside  diameter  of  152.4  mm. 
Immediately  after  the  extension  section  is  the  reduction  section.  It  is  0.61  m  long;  with 
it’s  inside  diameter  varying  linearly  from  152.4  mm  at  its  inlet  to  76.2  mm  at  the  exit.  A 
nozzle  with  a  76.2  mm  inlet  diameter  can  then  be  attached  directly  to  the  exit  of  the 
reduction  section.  In  the  present  experiments  a  converging  axisymmetric  nozzle  having 
an  exit  diameter  of  50.8mm  was  used.  The  nozzle  upstream  of  the  exit  was  designed 
using  a  fifth-order  polynomial  with  a  contraction  ratio  of  approximately  2.25.  The 
stagnation  pressure  and  temperature  were  held  constant  to  within  0.5%  of  its  nominal 
value  during  the  experiment. 
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The  jet  exhausts  into  an  anechoic  chamber  that  measures  5.2  m  wide,  5.8  m  long 
and  4  m  high.  The  walls,  ceiling,  and  floor  are  covered  with  high  temperature  305  mm 
deep  wedges  that  were  supplied  by  Eckel  Industries,  Inc.  The  wedges  are  composed  of 
#703  fiberglass  and  they  have  a  fiberglass  cloth  covering  followed  by  a  mesh  covering. 
There  are  sections  of  metal  grating  covering  the  wedges  on  the  floor  so  that  personnel 
movement  during  experimental  setup  is  possible.  The  chamber  has  provisions  for 
adequate  ventilation  for  easy  entrainment  of  the  ambient  air  into  the  chamber  during  jet 
operation.  The  exhaust  system  was  also  treated  to  minimize  any  noise  generated  by  the 
impinging  flow.  The  Inverse  Square  Law  was  used  to  determine  the  acoustic 
characteristics  of  the  chamber.  The  characterization  was  done  with  and  without  the  metal 
grating  in  place.  The  cut-off  frequency  of  the  chamber  is  about  500  Hz  with  the  grating 
in  place,  while  it  is  about  300  Hz  without  the  grating. 

2.1b.  Data  Acquisition  and  Automation 

There  are  a  number  of  different  conditions  that  need  to  be  monitored  during  an 
experiment  in  the  facility.  It  is  for  this  reason  that  an  integrated  data  acquisition  and 
control  system  was  implemented.  The  system  consists  of  three  PCs:  the  air  control  PC, 
the  DAQ  (data  acquisition)  PC,  and  the  burner  PC.  The  air  control  PC  uses  a  LabVIEW- 
based  program  to  control  the  air  pressure  entering  the  SUE  burner.  The  program  also 
monitors  the  static  pressure  and  stagnation  temperature  in  the  measurement  section  as 
well  as  the  ambient  pressure,  relative  humidity,  and  temperature  in  the  anechoic  chamber. 
All  of  this  data  is  displayed  for  easy  viewing  and  is  stored  for  future  reference. 

The  DAQ  PC  uses  a  MATLAB-based  program  to  acquire  acoustic  data.  The  PC 
has  3  National  Instruments  PCI-6110  high-speed  data  acquisition  boards  installed  in  it. 
Each  of  the  boards  can  simultaneously  sample  4  channels  and  each  channel  can  be 
sampled  at  a  rate  of  5  MS/s  with  12-bit  resolution.  Each  board  is  then  connected  to  its 
own  National  Instruments  BNC-21 10  connector  block.  The  MATLAB  program  can  be 
configured  to  sample  any  combination  of  the  12  available  channels  at  any  rate  up  to  the 
maximum.  The  program  also  sends  a  trigger  to  the  LabVIEW  program  on  the  air  control 
PC  so  that  all  of  the  relevant  conditions  are  logged. 

The  burner  PC  uses  a  Visual  BASIC-based  program,  supplied  by  Kaiser 
Marquardt,  to  control  the  SUE  burner.  It  monitors  the  various  pressures,  flow  rates,  and 
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temperatures  of  the  burner.  It  also  allows  for  the  control  of  the  following  solenoid  valves 
located  on  the  burner:  main  fuel,  igniter  fuel,  spark  plug,  and  nitrogen.  This  program  also 
controls  the  main  fuel  control  valve.  Due  to  the  importance  of  this  PC  it  is  isolated  from 
the  other  PCs  and  it  is  not  a  part  of  any  network. 

A  modem  acoustic  instrumentation  system  has  been  acquired.  This 
instrumentation  consists  of  ten  Bruel  &  Kjaer  model  4939  microphones  with  B&K  model 
2670  preconditioning  amplifiers  and  three  B&K  model  2690  Nexus  conditioning 
amplifiers.  The  signals  from  the  amplifiers  are  routed  through  a  breakout  box  that  leads 
to  the  three  BNC-2110  connector  blocks  mentioned  above.  Great  care  was  taken  to 
insure  that  all  of  the  cables  from  each  of  the  amplifiers  to  the  breakout  box  and  then  from 
the  breakout  box  to  the  connector  blocks  were  the  same  length.  The  microphones  were 
set  up  in  an  arc  that  had  a  radial  distance  of  50  diameters  from  the  nozzle  exit.  The  arc 
covered  the  polar  angle,  0,  range  from  90  to  150  degrees  relative  to  the  jet  inlet  axis. 

Each  of  the  microphones  had  a  relatively  flat  frequency  response  up  to  100  kHz 
and  was  subsequently  sampled  at  250  kHz.  The  dataset  for  each  microphone  contained 
409,600  samples  (1.6  sec).  This  allowed  for.  a  fast  Fourier  transform  (FFT)  of  4096 
points  over  100  subsets.  Averaging  the  results  for  the  100  subsets  reduced  the  random 
error  in  the  calculation  to  within  0.1%.  The  resulting  narrowband  spectrum  had  a  spectral 
resolution  of  61  Hz. 

The  overall  sound  pressure  level  (OASPL)  can  be  found  through  the  use  of  two 
methods.  The  first  method  involves  integrating  under  the  power  spectral  density  curve  to 
obtain  the  squared  pressure  fluctuation  and  then  applying  the  following  well  known 
formula: 

OASPL  =  101og^-  (1) 

Pref 

where  pref  is  20  pPa.  The  second  method  uses  the  pressure  data  that  is  obtained  from  the 
microphones  to  compute  the  root  mean  square  pressure.  This  value  is  then  squared  and 
then  Equation  1  can  then  be  used.  In  both  cases  the  OASPL  value  was  found  to  be  the 
same. 

A  number  of  corrections  must  be  applied  in  order  to  obtain  accurate  data.  The 
sound  pressure  level  (SPL)  at  each  frequency  needs  to  be  determined  through  the  use  of 
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the  raw  data.  The  corrections  for  the  actuator  response  as  well  as  the  free-field  response 
are  then  applied  at  each  frequency.  Lastly,  the  effect  of  atmospheric  absorption  at  each 
frequency  needs  to  be  determined  and  applied.  This  was  done  through  the  use  of 
formulas  provided  by  Blackstock14.  The  corrected  SPL  values  are  then  converted  back 
into  pressure  values  and  integration  is  then  performed  over  the  corrected  spectrum.  The 
resulting  squared  pressure  value  can  then  be  used  to  obtain  the  OASPL. 

To  illustrate  the  suppression  technique  on  a  full-scale  engine,  the  frequencies  are 
divided  by  a  factor  of  12  (assuming  the  full  scale  nozzle  diameter  =  0.6  m)  and  then  the 
frequency  spectra  are  converted  into  a  discrete  one-third-octave  spectrum.  In  order  to  do 
this,  one  must  first  determine  the  center  frequency  of  each  octave  band  and  then  its  lower 
and  upper  limits.  The  corrected  pressure  spectrum  is  then  integrated  over  these  limits  to 
determine  the  SPL  value  at  the  respective  center  frequency.  Then  the  entire  spectrum  is 
converted,  with  the  A-weighting  to  properly  reflect  subjective  judgments  of  the  noise  as 
commonly  used  in  the  literature. 


Figure  4.  The  Microjet  Arrangement  with  Respect  to  the  Primary  Nozzle.  Primary  Nozzle  Diameter 
=  50.8  mm.  Microjet  Nozzle  Diameter  =  400  pm. 


For  the  laboratory  scale  experiments  reported  here,  a  converging  axisymmetric 
nozzle  with  an  exit  diameter  of  50.8  mm  was  used.  The  micro-nozzles  were  made  of  400 
pm  stainless  steel  tubing.  The  underexpanded  microjets  impinge  on  the  shear  layer  at 
6mm  downstream  of  the  nozzle  exit.  The  angle  of  the  microjets  with  respect  to  the 
upstream  jet  centerline  was  60  degrees.  Typically  eight  microjets  are  used,  but 
experiments  were  also  carried  out  with  four  and  sixteen  microjets.  The  microjet 
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arrangement  is  shown  in  Figure  4.  The  stagnation  pressure  of  the  microjets  was  varied 
from  300  psia  to  700  psia.  At  500  psia,  the  fully  expanded  Mach  number  of  the  microjet 
is  about  2.91  and  the  corresponding  velocity  and  Reynolds  number  of  the  microjet  are 
6 18  m/s  and  about  8.7xl03,  respectively. 

The  jet  exhausted  into  a  quiet  surrounding  at  ambient  conditions.  The  stagnation 
temperature  of  the  primary  jet  was  kept  at  1033  K.  The  fully  expanded  jet  velocities  for 
the  Mj  =  1.38  is  760  m/s.  The  corresponding  Reynolds  number,  based  on  the  nozzle  exit 
diameter,  is  4.8  x  105. 

2.2.  Results  and  Discussion 

Figure  5  shows  shadowgraph  pictures  of  the  primary  jet  with  and  without  the 
microjet  injection  taken  using  a  high-speed  digital  camera  operating  at  1000  frames/sec 
with  an  exposure  time  of  28  (isec.  The  fully  expanded  Mach  number  of  the  primary  jet  is 
1.38.  As  expected,  a  bow  shock  in  front  of  the  microjet  appears  prominently  in  the 
picture  on  the  right.  A  faint  but  noticeable  trace  of  the  microjet  is  also  visible  in  the 
picture.  The  penetration  distance  of  the  microjet  is  known  to  be  primarily  a  function  of 
the  momentum  ratio15.  For  the  conditions  of  the  present  experiment  with  Mj  =  1.38,  the 
correlation14  gives  a  value  of  about  5mm,  which  is  close  (~  6mm)  to  the  measurement 
from  the  shadowgraph  visualization.  The  total  mass  flow  rate  of  the  microjets  at  500  psia 
is  about  1.1%  of  the  primary  jet  mass  flux. 


Figure  5.  Shadowgraphs  of  the  Jet  with  and  without  Microjet  Injection.  Left:  Normal  Jet;  Right: 

Air  Injection  at  500  psia 
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Figure  7.  Non-Dimensional  Centerline  Jet 
Temperature  with  and  without  Microjet 


Figure  6.  Infrared  Pictures  of  the  Mj  - 138  Jet. 

Top:  Normal  Jet;  Bottom:  Air  Injection  at  500 
psia 

Pictures  of  the  Mj  =  1.38  jet  taken  using  an  infrared  camera  with  and  without  the 
microjets  are  shown  in  Figure  6.  Due  to  the  underexpanded  nature  of  the  jet,  several 
shock  cells  are  seen  in  the  picture.  The  observed  effect  of  the  microjets  is  to  reduce  the 
shock-cell  length  and  the  corresponding  temperature  distribution  as  shown  in  Figure  7. 

Figure  8  shows  the  pressure  time  history  of  the  far-field  microphone  signal  at  the 
peak  radiation  angle  (0  =140°)  and  at  normal  to  the  jet  axis  for  Mj  =  1.38.  A  prominent 
characteristic  of  the  signal  at  the  peak  radiation  angle  is  the  random  occurrence  of  distinct 
bursts  of  strong  narrow  positive  pressure  transients  as  shown  in  Figure  9.  Such  a  signal  is 


Time  (sec)  Time  (sec) 


Figure  8.  Far-field  Microphone  Pressure  Time  Signals  for  a  Normal  Jet  (red)  with  Microjet  Injection 
(blue).  Left:  Peak  Radiation  Direction;  Right:  Normal  to  the  Jet  Axis. 


representative  of  a  crackling  jet,  first  recognized  and  described  by  Ffowcs  Williams  et 
al8.  Crackle  levels  tend  to  peak  at  near  the  eddy  Mach  angle  given  by  the  following 
equation5: 


sinor  =  —  ;M,  =^-^ 
M‘  &.  +  1 


where  Mj  and  a /  are  the  fully  expanded  Mach  number  and  speed  of  sound  of  the  primary 
jet,  respectively,  ax  is  the  ambient  speed  of  sound  and  a  is  measured  clockwise  from  the 
jet  axis.  The  estimated  peak  radiation  angle,  0,  using  Equation  2  is  139°,  which  is  in 
agreement  with  the  present  results  within  the  experimental  uncertainty.  Also,  included  in 
Figure  7a  is  the  signal  obtained  with  the  microjet  injection;  the  lack  of  the  high  amplitude 
positive  peaks  with  injection  is  clearly  depicted,  suggesting  crackle  suppression.  Figure 
8b  shows  the  signals  taken  at  90°  with  and  without  the  microjet  injection.  In  this 
particular  trace,  the  instant  at  which  the  microjet  injection  turned  on  is  also  noted.  The 


Figure  9.  Typical  Far-field  Microphone  Crackle  Signal. 


relatively  low  crackle  levels  found  in  the  signal  at  this  angle  of  observation  are  in 
conformity  with  the  previous  observations8,10.  The  effect  of  microjets  on  the  signal  at 
normal  to  the  jet  axis  is  mostly  reflected  in  the  reductions  of  the  overall  signal  strength. 
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Ffowcs  Williams  et  al8  have  suggested  that  the  skewness  factor  of  the  recorded 
signal  is  an  effective  direct  measure  of  crackle.  Keeping  this  in  mind,  the  skewness  of  the 
far  field  microphone  signal  at  different  angles  with  and  without  microjet  injection  is 
shown  in  Figure  10.  The  skewness  for  jets  peaks  near  the  peak  radiation  angle  when  the 
eddy  Mach  number,  Me,  is  greater  than  1.2.  Away  from  the  peak  radiation  angle,  the 
skewness  drops  rapidly,  an  indicator  of  the  absence  of  crackle.  Indeed,  the  signal  at  90°, 
shown  in  Figure  7b,  confirms  this  finding.  The  skewness  values,  and  their  dependence 
on  the  angular  position  found  in  the  present  experiment,  are  in  general  agreement  with 
those  of  Ffowcs  Williams  et  al8.  In  support  of  the  earlier  observation,  with  respect  to  the 
pressure  time  signals,  the  crackle  suppression  in  the  aft  quadrant  of  the  jet  exhaust  is 
quite  evident  by  skewness  reductions.  The  energy  contained  in  the  spikes  responsible  for 
crackle  is  about  30%  of  the  total  energy7’10,  as  such  it  is  expected  that  the  OASPL  in  the 
peak  radiation  angle  direction  be  reduced  by  at  least  1 .5  dB. 


Figure  10.  The  Directional  Distribution  of  the  Skewness  Factor  with 
and  without  Microjet  Injection.  Mj  =  1.38. 

Figure  1 1  shows  the  variation  of  the  OASPL  with  0  for  the  normal  jet  and  with 
microjet  injection.  Significant  reductions  in  the  aft  quadrant  are  observed.  At  the  peak 
radiation  angle,  the  reduction  amounts  to  about  4.5  dB.  Since,  crackle  is  known  to 
contribute  only  about  1.5  dB,  the  additional  reductions  in  the  OASPL  can  be  attributed  to 
suppression  of  part  of  the  unsteady  flow  generated  by  the  large  structures,  commonly 
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refereed  to  as  large-scale  mixing  noise.  Measurements  taken  at  Mj  =  0.9  jet,  discussed 
later  suggests  that  the  large-scale  mixing  noise  is  also  suppressed  to  some  extent  by  the 
microjet  injection. 


Figure  1 1 .  Far-field  Directivity  of  Mj  =  1.38  with  and 
without  Microjet  Injection.  Nozzle  Temperature  Ratio  =  3. 


Spectral  measures  of  the  microphone  signals  show  the  frequencies  affected  by  the 
microjet  injection  and  thereby  lead  to  some  observations  about  the  influenced  source 
types.  The  spectrum  of  a  sharp-edged  spike  shown  in  figure  8  is  flat  and  evenly 
distributed  over  a  wide  range  of  frequencies.  Hence,  the  crackle  suppression  should  lead 
to  reductions  in  a  wide  band  of  frequencies  in  the  spectrum.  Figure  12  shows  the  narrow 
band  (Af  =  61  Hz)  frequency  spectra  corresponding  to  the  far-field  microphone  signals 
shown  in  Figure  8.  Although,  the  computed  spectra  extends  up  to  100  kHz,  for 
compactness  only  the  data  up  to  60  kHz  is  shown  here.  Screech  tones  that  are  commonly 
present  in  underexpanded  are  clearly  seen  in  the  spectra  with  its  fundamental  tone  being 
at  3662  Hz.  The  corresponding  Strouhal  number  (St  =  fD/U)  of  0.26  is  in  agreement 
with  previous  measurements  reported  in  the  literature16.  The  suppression  of  screech  tones 
due  to  microjet  injection  is  clearly  evident.  But  more  importantly,  at  peak  radiation  angle, 
the  SPL  reductions  are  observed  in  the  entire  spectrum  (Figure  12a).  Except  for  the 
absence  of  screech  tones,  the  spectrum  in  Figure  12b,  corresponding  to  0  =  90°,  shows 
little  change  with  microjet  injection.  Since,  the  conventional  fine  and  intermediate  scale 
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Figure  12.  Narrowband  Frequency  Spectra  for  the  Mj  =  1 .38  Jet.  Left:  Peak  Radiation  Angle;  Right: 

Normal  to  the  Jet  Axis. 


Figure  13.  A-weighted  Spectra  Extrapolated  to  the  Full  Scale  Corresponding  to  Figure  1 1 . 


chaotic  turbulence  contributes  mostly  to  the  noise  at  this  angle,  it  is  suggested  that  the 
microjets  have  little  influence  on  this  component  of  the  noise.  Additionally,  the 
broadband  shock-associated  noise  that  is  dominant  in  the  forward  quadrant  and  at  90o  is 
also  minimally  affected.  From  these  observations,  it  is  evident  that  the  microjet  injection 
influences  the  crackle  and  the  large-scale  mixing  noise.  The  source  of  low  amplitude 
undulations  seen  in  the  spectra  is  currently  being  investigated,  but  they  seem  to  have  little 
influence  on  the  OASPL  based  on  the  agreement  of  the  present  data  with  that  of  Tanna  et 
al17. 

Using  the  scaling  factor  of  12  to  account  for  the  full-scale  nozzle  and  the  human 
response  to  noise,  the  narrowband  spectra  are  recalculated  to  yield  A-weighted  spectra. 

16 


Shown  in  this  fashion  in  Figure  13,  at  the  peak  radiation  angle,  the  SPL  reductions  are 
quite  significant  (almost  10  dBA)  except  at  very  low  frequencies  (less  than  500  Hz).  The 
corresponding  reduction  in  the  OASPL  is  about  6  dBA.  The  reduction  of  1.5  dBA  in  the 
OASPL  at  0  =  90°  may  mainly  be  due  to  the  screech  elimination. 
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3.  ACTIVE  NOISE  CONTROL  ARCHITECTURE 

3.1.  Introduction 

The  control  system  addressed  in  this  project  addresses  the  active  quieting  of  jet 
noise  as  shown  below  in  Figure  14.  The  developments  of  the  program  documented  herein 
were  designed  to  address  the  emanated  noise  by  introducing  a  modulated  flow  of  cold  air 
into  the  nozzle  just  prior  to  exit  in  such  a  fashion  as  to  create  canceling  noise  in  the  far 
field.  Sensors  of  pressure  would  be  used  inside  the  nozzle  to  sense  incipient  fluctuations, 
and  these  measured  pressure  variations  would  be  used,  with  appropriate  adaptive 
filtering,  to  drive  the  modulation  actuators  to  achieve  attenuation  in  the  far  field.  If 
necessary,  a  towed  array  of  pressure  sensors  (microphones)  would  be  used  to  help 
characterize  and  optimize  the  quieting. 


The  program  had  two  main  thrusts:  development  and  demonstration  of  a  flow- 
modulation  actuator,  and  implementation  and  demonstration  of  an  adaptive  noise- 
cancellation  algorithm.  Actuator  development  was  based  on  PSI’s  (Planning  Systems 
Inc.)  previously  proposed  concept  of  a  fast  piezo-operated  valve.  Originally,  four 
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actuators  were  envisioned.  However,  when  a  more  thorough  investigation  of  piezo 
elements  was  undertaken,  it  was  found  that  only  high-voltage,  high-capacitance  elements 
would  be  capable  of  the  continuous-duty  displacements  at  1000  Hz.  The  amplifiers 
necessary  to  drive  such  elements  were  much  more  costly  than  estimated.  Further,  the 
difficulty  of  valve  element  implementation  forced  a  more  step-wise  development.  Two 
piezo  elements  were  purchased,  with  one  high-current  amplifier.  Two  separate  versions 
of  the  actuator  were  implemented,  as  will  be  discussed  below.  The  first  of  these 
encountered  two  main  problems  which  hampered  its  results: 

•  The  valve  elements  were  not  successfully  machined  by  the  vendor,  causing  delays 
in  implementation,  excessive  leakage  and  friction  and  sub-standard  friction  and 
closing  performance. 

•  The  concept  of  the  actuator  caused  excessive  energy  loss,  due  to  the  small  hole 
sizes,  prior  to  injection  into  the  hot  jet. 

However,  the  first  actuator  did  demonstrate  that  the  piezo  element  could  be  successfully 
integrated  to  drive  the  element  at  full  deflection  up  to  1000  Hz.  A  revised  actuator  was 
manufactured  that  increased  the  opening  area  by  a  factor  of  three,  and  largely  resolved 
the  tolerancing  and  hole  mismatching  problems.  In  addition,  a  small  plenum  was  added 
after  the  valve  element  to  reduce  the  energy  loss  prior  to  injection  into  the  main  jet.  Both 
low-frequency  flow  modulation  and  energy  modulation  were  as  designed.  However,  use 
of  the  secondary  plenum  precluded  operation  at  acoustic  frequencies,  because  the  plenum 
smoothed  out  and  attenuated  high-frequency  flow  variations.  A  third  actuator  was  then 
proposed  with  direct,  large-hole  modulation  via  a  levered  piezo  device.  However, 
funding  was  not  available  to  implement  it.  The  adaptive  algorithm  was  developed  and 
implemented  essentially  as  proposed.  However,  since  no  effective  actuator  was  available 
for  testing  with  the  miniature  jet,  demonstration  in  PSI’s  laboratory  of  the  algorithm 
function  was  undertaken.  Approximately  10  dB  of  attenuation  of  a  colored  noise 
between  100  and  200  Hz  was  demonstrated. 

The  algorithm  has  the  following  capabilities: 

•  Multiple-sensor,  multiple-actuator  control. 

•  Feedback,  feedforward,  or  mixed  operation. 

•  Infinite-horizon  control  optimization. 
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•  No  prior  or  off-line  modeling  required. 

The  algorithm  was  implemented  in  two  major  components.  The  first,  which  embodied 
all  learning  and  adaptation,  was  implemented  in  Matlab  on  the  host  PC.  The  second, 
which  embodied  the  real-time  control  and  filtering  computations,  was  implemented  C 
language  on  a  C60  digital  signal  processor  residing  on  a  plug-in  card. 

3.2  Actuator  Description  and  Demonstration  Results 

3.2a.  Initial  Actuator  Design 

The  actuator  design  was  based  on  a  moving  valve  concept,  with  the  prime  mover 
a  piezo  element.  Previous  electromechanical  designs  have  been  limited  by  the  need  to 
drive  immense  currents  to  achieve  the  extreme  accelerations  required  by  high-frequency 
operation.  Such  high  currents  required  complicated  cooling  systems.  Piezo  elements 
have  the  advantage  of  the  innate  ability  to  achieve  extremely  large  accelerations,  albeit 
over  relatively  small  stroke.  Our  approach  to  the  valve  design  was  initially  fairly 
conventional,  except  for  the  tiny  hole  sizes  required  by  the  small  stroke  of  the  piezo 
element. 

Design  sizing  of  the  valve  opening  area  followed  the  analysis  by  Blondel  and 
Elliott21.  Design  of  the  actuator  centers  around  a  pressurized  plenum  chamber  separated 
from  the  source  output  by  a  valve.  Initially,  our  design  was  aimed  at  achieving  a  specific 
mass  flow  rate  modulation  of  air,  without  regard  to  the  energy  of  the  flow,  because  we 
were  not  aware  of  the  importance  of  this  aspect  of  the  actuator  to  its  modulation 
capabilities  of  the  jet-engine  exhaust.  Figure  15  illustrates  the  design  concept.  The  valve 
was  to  be  of  sufficient  size  to  allow  0.014  kg/sec  of  air  to  pass  through  the  valve  opening 
with  the  plenum  chamber  pressure  set  at  620  KPa. 
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Figurel5.  Electropneumatic  Transducer  Concept 


From  reference  21,  using  air  as  the  fluid  medium,  the  equation  needed  to  size  the  orifice 
for  a  given  set  of  conditions  is: 


Si En  a 
A  A 


where: 

p2  =  pressure  at  station  2 

vy  =  particle  velocity  at  station  2 

c0  =  speed  of  sound  in  the  medium 

ppt  =  pressure  in  the  plenum  chamber 

4  =  time  varying  area  or  valve  cross  sectional  area 

4  =  cross  sectional  area  of  the  output  section 

Air  Valve  Size 

The  cross  sectional  area  needed  for  the  valve  was  determined  as  follows.  We  know  for  a 
control  volume: 

P  =  4V2P2  (4) 


where  p  is  the  mass  flow  rate,  4  is  the  cross  sectional  area  of  the  output  section,  v2  is 
the  particle  velocity  at  Station  2,  and  p2  is  the  fluid  density  at  Station  2. 

So  combining  the  two  equations  we  have: 
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•  ,  fsY  coPpl.) 

p  =  A2v2p2=  -  — M  P2 

A  oy  p2  ) 


From  Boyle’s  Law: 

Pi  =  P2^2 

where  p2  is  the  pressure  at  Station  2,  R  is  the  universal  gas  constant,  and  T2  is  the 

absolute  temperature  downstream  of  the  valve.  By  substituting  —  =  RT2  into  our 

Pi 

equation  for  p  above  and  solving  for  Ax  we  get: 

,  .  rt2 
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For  our  initial  test  configuration  to  be  tested  at  the  PSI  Melbourne  laboratoiy  the 
following  values  were  used  to  calculate  Ax : 

p  =  0.014  kg!  sec  -  , 

R  =  287  Joules /kg°K 
T2=-40°C 
c0  =  340  meters!  sec 
ppl  =621  KPa 

From  these  values  we  calculated  Ax  to  be: 

Ax  =  7.75x1 0"6  meters 2  or  0.0 1 1  inches 2 

From  the  cross  sectional  area  we  determine  the  diameter  of  the  valve  to  be  0. 12 1  inches. 
Actuator  Drive 

Because  of  the  high  accelerations  associated  with  opening  and  closing  a  valve  at  up 
to 1000  Hz  weight  of  the  moving  mass  and  length  of  stroke  had  to  be  minimized.  We 
could  not  purchase  or  develop  an  actuator  capable  of  opening  and  closing  a  0.121  inch 
diameter  valve  at  1000  Hz.  Other  key  issues  influencing  the  choice  of  a  suitable  actuator 


•  Small  size 

•  Actuator  power  requirements 

•  Survival  near  a  jet  nozzle  (high  Temperatures) 

•  Self  generated  heat/survivability 

•  -  High  accuracy 

•  Linear  dynamic  response 


22 


A  piezo-actuator  was  chosen  as  the  only  suitable  means  to  drive  the  air  vane.  The  largest 
piezo-actuator  available  capable  Of  continuous  full  stroke  operation  at  1000  Hz  is  one 
with  0.007  inch  travel. 

Air  Valve 

A  series  of  0.007  inch  diameter  holes  were  used  to  realize  the  cross  sectional  area  of 
0.0 1 1  inches2  while  the  valve  is  open.  For  an  open  area  of  0.0 1 1  inches2  we  needed  a 
vane  with  approximately  292  0.007”  diameter  through  holes.  Two  concentrically  located 
thin  walled  stainless  steel  tubes  make  up  the  air  vane.  To  achieve  symmetry  the  vanes 
were  manufactured  with  12  rows  of  24  holes  or  288  holes  total.  Figure  16  shows  the 
concentric  tube  air  vane.  The  inner  vane  is  stationary  while  the  piezo-actuator  drives  the 
outer  vane. 


Figure  16.  Concentric  Tube  Air  Vane 

Flexures 

Logarithmic  spiral  linear  flexure  bearings  were  used  to  position  and  guide  the  vane.  The 
flexures  were  manufactured  out  of  Beryllium  Copper  0.010  inches  thick.  A  flexure  was 
attached  to  each  flange  of  the  moving  vane  to  insure  linearity  of  movement.  The  flexures 
and  positioning  hardware  are  shown  in  17  attached  to  the  vanes. 
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Figure  17.  Flexures  and  Positioning  Hardware 


Actuator/Air  Vane  Assemble 


The  assembled  actuator  is  shown  in  Figure  18.  The  piezo-actuator  is  positioned  in-line 
axially  to  the  air  vanes.  A  moment  de-coupling  flexure  is  mounted  between  the  actuator 
and  the  vane  to  eliminate  any  off  axis  moments  introduced  to  the  piezo-actuator  during 
operation.  Off-axis  moments  can  crack  the  ceramics  on  the  actuator. 

Plenum  Chamber 

Figure  19  shows  the  plenum  chamber  with  the  actuator  and  associated  positioning  and 
spacing  brackets  assembled.  With  the  actuator  assembly  mounted  inside  the  plenum 
chamber  a  micrometer  is  placed  in  series  with  the  inner  vane  and  piezo-actuator  to  allow 
adjustment  of  the  vane  to  a  fully  closed  position  prior  to  modulation.  The  outer  vane  is 
fixed  at  a  known  location.  With  a  constant  pressure  air  source  the  plenum  chamber 
maintains  the  quasi-static  positive  air  pressure.  The  end-caps  are  used  to  seal  the 
pressure  vessel  and  operate  as  a  base  to  position  the  interior  components.  The  plenum 
chamber  is  fabricated  out  of  Plexiglas  so  we  can  monitor  the  moving  components.  The 
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end-caps  are  Aluminum.  The  figure  shows  the  plenum  chamber  with  the  actuator  and 
associated  positioning  and  spacing  brackets  assembled. 


Air  Vane  to 


Initial  Actuator  Test  Setup 


The  actuator  test  setup  is  shown  in  Figure  20.  Two  Honeywell  4000  pressure  transducers 
were  used  to  measure  plenum  chamber  and  downstream  pressure.  A  Piezotronics 
pressure  transducer  with  a  1  KHz  sample  rate  was  also  used  downstream  of  the  valve.  A 
Flowstream  Type  4  OFS  flow  sensor  was  used  to  measure  the  modulated  flow  rate.  Test 
results  using  these  sensors  in  the  test  setup  shown  are  presented  in  the  following  section. 
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Figure  19.  Plenum  Chamber  with  Actuator  and  Air  Vanes 


Piezotronlcs  Pressure  Sensor 
(1  KHz  Sample  Rate) 


Honeywel  4800PC  Series 
Pressure  Transducers 


Rowstream  Type  4  OFS 
.Series  How  Sensor 


r  OD  *  12*  Long  Pressure  Vessel 


Air  Pressure  Source 


Figure  20.  Actuator  Test  Setup 
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3 ,2b.  Initial  Actuator  T est  Results 

The  initial  actuator  was  tested  during  a  2-day  test  visit  at  Dr.  Krothapalli’s  FSU 
facility.  The  purpose  of  the  testing  was  to  measure  the  performance  of  the  actuator  when 
mounted  on  the  jet  nozzle.  The  actuator  was  attached  to  the  jet  nozzle.  Cables  of  several 
sensors,  integral  to  the  actuator,  and  the  piezo  element  provide  interface  from  the  actuator 
to  the  data  acquisition  system.  Also,  several  other  FSU-provided  sensors,  including 
Kulite  pressure  sensor  (mounted  at  the  jet  nozzle),  and  B&K  microphones  positioned  at 
the  near  and  far-field,  were  monitored  by  the  data  acquisition  system.  During  the  testing, 
a  flow  rate  sensor  was  added  upstream  to  the  air  source  of  the  actuator.  This  sensor’s 
output  was  also  monitored  by  the  data  acquisition  system.  The  Table  1  lists  all  the  test 
scenarios  that  were  performed. 
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Table  1  Description  of  Test  Cases  and  Results  for  FSU  Test  Visit 


Test  Case 
# 

Description  of  Test  Case 

Significant  Results 

initial  actuator  checkout 

1  thru  5 
(day  1) 

Tests  performed  to  determine 
displacement  of  actuator’s  vane. 

20, 200, 1000  Hz  sine,  full  voltage  drive 
to  piezo,  measuring  actual  applied 
voltage  and  position  displacment  of 
actuator’s  vane. 

Only  getting  ±1 .4  mils  of  motion, 
when  expecting  ±3.5  mils.  Also, 
motion  is  highly  distorted.  We  blame 
binding  (rubbing)  of  inner  and  outer 
vane  for  both  problems. 

Refer  to  Figure  21  for  typical  vane 
displacement  vs  time  plot. 

_ '  _ sound  pressure  testing _ 

6  and  1 6  These  tests  were  performed  to  determine  No  difference  in  sound  pressure  (at  the 
sound  pressure  due  to  actuator  flow.  excitation  tone)  was  measured  with  and 

(with  no  flow  in  the  jet  nozzle)  without  air  flow  thru  the  actuator.  This 

is  indicative  of  the  modulated  flow  rate 
Monitored  Kulite  pressure  and  ‘nearby’  being  much  less  than  expected, 
outside  microphone  (in  addition  to  piezo 

voltage  and  vane  displacement).  Under  Refer  to  Figure  22  for  plots  of  ambient 
following  conditions:  noise,  20  Hz  piezo  drive  -  no  actuator 

•  ambient -no  flow,  no  actuator  drive  flow,  and  20  Hz  piezo  drive  -  with  40 

•  no  flow  -  actuator  set  a  mid  stroke  Ps*  *n  actuator  chamber. 

•  no  flow  -  actuator  driven  at  20, 200 
and  1000  Hz  at  full  level 

•  40  psig  flow  -  actuator  not  driven 

•  40  psig  flow  -  actuator  driven  at  20, 

200  and  1000  Hz  at  full  level 

1 7  thru  20  These  tests  were  performed  to  determine  Similar  results  as  previous  tests.  No 
the  sound  pressure  variation  due  to  distinguishable  sound  pressure 

actuator  flow  while  there  is  (cold)  flow  variation  due  to  actuator  flow  output 
in  the  jet.  was  seen  in  near-field  microphone 

response. 

Also,  added  a  data  acquisition  channel  to 
monitor  an  upstream  (to  the  actuator) 
flow  rate  sensor  (in  addition  to  other 
sensors). 

•  Re-ran  some  previous  tests. 

•  85  psig  flow  in  actuator. 

•  44  psig  flow  thru  jet  (cold) 

21  thru  23  Ran  several  tests  to  determine  if  Measured  86  dB  SPL  at  nearby 

(day  2)  wrapping  ‘Cerablanket’  material  around  microphone  with  and  without 

actuator  housing  would  lessen  acoustic  ‘Cerablanket’  material.  Conclusions: 
transmission  of  noise  due  to  piezo  material  had  no  effect  on  reducing  this 

motion.  _  undesirable  noise. 
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Test  Case 
# 

Description  of  T est  Case 

Significant  Results 

24  thru  26 

Ran  several  tests  to  determine  effect  of 
increase  air  pressure  in  actuator  on  flow 
rate  modulation. 

•  5  Hz,  full  level  drive  to  piezo 

•  tested  40, 72  and  1 1 2  psig  (could  not 
get  higher  pressure  due  to  pressure 
drop  between  air  tank  and  actuator, 
and  we’re  concerned  about  pressure 
limitation  of  actuator) 

The  mean  flow  rate  increase  from  180 
to  480  slpm,  while  the  modulated  flow 
rate  at  5  Hz  increase  by  factor  5.5 
(when  increasing  40  to  122  psig 
actuator  inlet  pressure). 

Refer  to  Figure  23  for  PSD  of  flow  rate 
sensor. 

effect  of  increasing  actuator  plenum  pressure 
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Added  far-field  microphone  to  test  set¬ 
up.  Recorded  this  and  other  sensors  with 

•  44  psi  cold  flow  in  jet 

•  122  psi  in  actuator 

Also,  attempted  to  run  ‘hot’ flow  jet,  but 
were  unable  to  light  the  burners. 

We  were  able  to  show  fairly  good 
correlation  between  Kulite  pressure 
and  the  ‘far-field’  microphone,  which 
is  gives  us  hope  in  being  able  to  use  the 
Kulite  sensor  as  feedback  in  the  noise 
cancellation  algorithm.  Refer  to  Figure 
24  for  correlation  result. 

actuator  flow  rate  modulation 

Ran  several  tests  to  determine  static  flow  I  plenum  pressure  DC  &  modulated 


of  actuator  vs.  actuator  chamber 
pressure. 

flow 

40  psig 

72  psig 

122  psig 

165  &  20  SLPM 
270  &  30  SLPM 
450  &  45  SLPM 

actuator  stagnation  pressure 

Tests  were  run  to  determine  the 

plenum  pressure 

outlet  pressure 

stagnation  pressure  at  the  outlet  of  the 
actuator  vs.  actuator  chamber  pressure. 

40  psig 

72  psig 

122  psig 

0.35  psig 
0.875  psig 
2.55  psig 

These  tests  verified  that  actuator  vane 
design  would  not  provide  outlet 
pressure  needed  to  modulate  far-field 
sound. 

Far  Field  PSD  (p$r/H2) 
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flowrate  measurement  of  runs  24,  25 ,26 


Figure  23.  PSD  of  Flow  Rate  Sensor  for  Varying  Actuator  Inlet  Pressure 


F igure  24.  Correlation  of  Kulite  Pressure  T ransducer  to  F ar-F  ield  Microphone 


To  address  the  low  flow  rate  modulation  issue,  we  formulated  a  plan  that  would 
increase  it  by  5x.  This  plan  included  using  our  improved  vane  design  (with  more  holes, 
and  fabricated  by  a  higher  precision  machine  shop).  From  several  discussions  with  Dr. 
Krothapalli  during  the  final  day  of  the  test  visit,  it  was  determined  that  the  approach  of 
the  initial  actuator  design  was  misguided.  The  actuator  requires  larger  holes  (instead  of 
many  small  ones)  that  are  closer  to  the  flow  stream  of  the  jet.  The  next  section  describes 
modifications  made  to  the  initial  design  to  incorporate  this  revised  approach. 

3.2c.  Revised  Actuator  Design 

Results  of  testing  performed  on  the  Initial  Actuator  proved  the  cross  sectional  area 
of  the  valve  was  too  small.  Irregularities  in  the  fabrication  work  contributed  to  the  flow 
rate  through  the  vane  being  lower  than  expected.  The  mass  flow  rate  was  approximately 
10%  of  what  was  needed  in  a  full  open  position.  The  new  valve  design  more  than  tripled 
the  total  cross  sectional  area.  The  revised  vane  had  912  0.007  inch  holes  in  19  rows  of  48 
holes.  The  total  open  area  increased  from  0.0 1 1  inches 2  to  0.035  inches 2 .  In  addition  to 
increasing  the  open  area,  stringent  dimensioning  and  tolerancing  methods  were  used  to 
insure  the  new  vane  was  fabricated  as  designed.  In  addition  to  changing  the  open  flow 
area  of  the  vane,  we  integrated  a  Stainless  Steel  plenum  chamber  to  house  the  actuator. 
The  plenum  chamber  was  designed  to  interface  to  the  ports  on  the  Florida  State 
University  anechoic  chamber  jet  nozzle.  A  thermal  barrier  plate  was  included  between 
the  plenum  chamber  and  the  nozzle  because  of  the  extremely  high  temperatures  on  the 
surface  of  the  nozzle  seen  during  operation.  A  venturi  effect  orifice  plate  was  added 
between  the  jet  nozzle  and  the  thermal  barrier  plate  to  magnify  the  modulation  effect  of 
the  vane.  By  introducing  this  plate  between  the  valve  and  the  source  output  more  energy 
is  transferred  to  the  flow  stream.  Figure  25  illustrates  the  new  plenum  chamber,  thermal 
barrier  plate  and  orifice  plate.  Figure  26  shows  a  cross  section  of  the  interface  plates 
installed  on  the  FSU  jet  nozzle. 
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Figure  26.  Actuator  Interface  to  Jet  Nozzle 
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3 .2d.  Revised  Actuator  Test  Results 

The  final  actuator  was  tested  in  the  lab  facilities  of  PSI's  Melbourne,  FL  office. 
Unfortunately,  several  last-minute  component  failures  occurred  that  prevented  testing  at 
the  FSU  facility.  As  mentioned  previously,  the  initial  actuator  design  was  modified  to 
include  an  intermediate  chamber  to  increase  the  outlet  pressure.  This  intermediate 
chamber  and  the  outlet  pressure  were  measured  for  various  outlet  hole  sizes,  actuator 
input  pressures  and  actuator  drive  frequencies.  Figure  27  shows  the  testing  configuration. 
Table  2  below  shows  the  results  of  various  test  cases.  The  inlet  pressure  was  measured 
with  a  pressure  gauge,  the  outlet  pressure  was  measured  using  a  pitot  tube  positioned  at 
the  exit  of  the  actuator’s  intermediate  chamber,  and  the  flow  rate  was  measured  using  an 
electronic  flow  meter  (Univeral  Flow  Monitor’s  model  number  OFS-MEBS700SLMN-6- 
A-X7A-D3). 


Intermediate 

Chamber 
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Table  2  Final  Actuator  Test  Cases  and  Results 


Test  Case 
# 

Description  of  Test  Case 

Results 

1 

Hole  restriction:  none 

Actuator  inlet  pressure:  40  psig 

Drive  frequency:  DC 

Actuator  Outlet  Pressure: 

Actuator  Open:  8.0  psi 

Actuator  Closed:  1.2  psi 

Actuator  Flow  Rate: 

Actuator  Open:  510  slpm 

Actuator  Closed:  260  slpm 

2 

Hole  restriction:  1/8”  diameter 

Actuator  inlet  pressure:  40  psig 

Drive  frequency:  DC 

Actuator  Outlet  Pressure 

Actuator  Open:  34.8  psi 

Actuator  Closed:  18.4  psi 

Actuator  Flow  Rate: 

Actuator  Open:  350  slpm 

Actuator  Closed:  230  slpm 

3 

Hole  restriction:  1/8”  diameter 

Actuator  inlet  pressure:  40  psig 

Drive  frequency: 

200  Hz 

100  Hz 
'20  Hz 

10  Hz 

Actuator  Outlet  Pressure: 

rms  level  measure  at  drive  frequency 

1.5  psi 

2.6  psi 

3.0  psi 

2.5  psi 

These  results  show  that  adding  the  intermediate  chamber  does  produce  the  desired  result 
of  increasing  the  stagnation  (or  outlet)  pressure.  Further  testing  was  disrupted  due  to 
failure  (cause  unknown)  of  the  piezo  actuators  and  amplifier. 


3.3.  Control  Algorithm  Description  and  Results 

The  control  algorithm  takes  the  measured  noise  signals  (from  sensors  at  the 
engine  nozzle  exit,  in  the  near  field,  or  in  the  far-field)  and  transforms  them  dynamically 
to  appropriate  drive  signals  for  the  actuators.  Due  to  the  unpredictable  and  complex 
nature  of  the  relationship  between  the  actuator  signals  and  the  measurements,  the  control 
algorithm  developed  and  demonstrated  by  PSI  is  adaptive  in  nature:  it  constantly 
measures  and  updates  its  model  of  the  system  to  be  controlled  and  changes  the  control 
filtering  accordingly. 

The  control  algorithm  is  based  on  an  architecture  used  by  PSI  for  several  years. 
The  system  identification  component  was  previously  developed  on  an  SBIR  Phase  I 
project  for  the  Air  Force  Research  Lab.  The  splitting  of  the  implementation  of  the 
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algorithm  into  a  PC-host-Matlab  component  for  heavy  computations  and  executive 
control,  and  a  DSP-target-C  component  for  real-time  computation  is  new  to  this  project. 
We  first  describe  the  control  system  at  the  top  level,  then  proceed  to  more  detailed 
information  on  the  theory  and  its  implementation.  A  laboratory  demonstration  of  the 
software  and  hardware  was  conducted,  and  is  described  last,  in  Section  0 
3.3a.  Control  Algorithm  Top  Level  Description 

The  noise  control  algorithm  is  a  general  adaptive  control  algorithm  that  is  capable 
of  both  feedforward  and  feedback  control  in  MIMO  systems.  It  addresses  the  control  of  a 
system  described  as  shown  below  in  Figure  28.  The  system  is  assumed  to  be  linear  time 
invariant,  and  all  quantities  are  assumed  to  be  vectors  in  discrete  time  (i.e.,  as  seen  by  the 
control  computer).  In  the  figure,  the  blue  boxes  represent  the  external  world,  and  the  red, 
the  controller  to  be  designed. 


Figure  28.  Description  of  Systems  Controllable  via  the  Adaptive  Algorithm 
System  disturbances  include  all  noises  (acoustic  or  otherwise)  that  drive  the 
system  and  have  significant  dynamics  (coloring)  as  seen  in  the  measurements.  Actuator 
commands  are  the  drive  outputs  sent  by  the  controllers.  The  controlled  system  includes 
the  dynamics  of  the  actuators  and  the  sensors  in  addition  to  those  of  the  jet  engine  or 
other  system  to  be  controlled. 

The  performance  errors  make  up  the  vector  to  be  minimized  by  control.  For  the 
far-field  noise  problem,  they  would  ideally  be  readings  of  the  far-field  instantaneous 
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noise  pressure.  However,  for  the  adaptive  control  system  they  must  be  accessible  signals, 
such  as  measurements  of  perturbations  near  the  exit  of  the  duct,  or  from  a  towed 
microphone  array.  In  the  control  setup,  they  are  assumed  to  be  a  known,  linear 
proportion  of  one  or  more  of  the  sensor  signals  (probably  microphones,  in  the  current 
instance),  prior  to  corruption  by  additive  white  noise  in  the  system. 

Note  that  this  setup  is  quite  general  and  includes  the  possibility  of  high-gain 
feedback  and  feedforward  control.  Feedforward  control  emerges  naturally  when  some  of 
the  sensor  readings  are  unaffected  by  the  actuator  commands,  in  which  case  they  can  be 
used  to  help  deduce  the  system  disturbances.  Note  also  that  by  judicious  choice  of  the 
setup  we  can  feed  through  some  of  the  actuator  signals  to  the  performance  errors — 
effectively  causing  a  tradeoff  to  be  introduced  between  the  level  of  actuation  effort  and 
the  other  performance  errors  of  the  system.  This  has  been  implemented  in  the  control 
system.  This  setup  is  precisely  that  of  the  classical  Linear  Quadratic  Gaussian  (LQG) 
problem.  In  fact,  we  have  chosen  to  use  LQG  to  design  the  controller  adaptively  in 
response  to  changes  in  the  identified  plant,  as  will  be  seen  in  greater  detail  below.  Our 
adaptive  control  algorithm  (the  red  block  of  the  previous  figure)  consists  of  the  principal 
blocks  depicted  in  Figure  29  below. 


Fast  Time  Computations 
(DSPS'QQOffz) 


probe  frequency 
emphasis  gains 


observer  | 

Stow  Time  Computations  prediction 

(Host/MatM)  error  history 

Figure  29.  Principal  Blocks  of  the  Adaptive  Noise  Control  Algorithm 
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The  computations  are  split  into  two  major  sections:  a  fast  set  of  computations, 
which  must  be  completed  at  the  full  sample  rate  of  the  system,  and  a  slower  set,  which 
contains  all  of  the  most  burdensome  calculations  associated  with  learning  and  adaptation. 
The  fast  set  is  programmed  in  C  code  and  runs  on  the  target  Texas  Instruments  C60 
processor,  whereas  the  slow  set  is  programmed  in  Matlab,  and  runs  on  the  host.  The 
setup  and  implementation  of  the  hardware  and  communications  is  discussed  further 
below. 

In  an  ideal  system  (all  dynamics  known  in  advance),  no  adaptation  would  be 
necessary,  and  the  whole  system  would  consist  just  of  the  controller  itself.  The  rest  of  the 
system  is  designed  to  ensure  that  the  correct  gains  and  parameters  are  used  in  the 
controller.  We  have  implemented  an  observer-based  controller.  Essentially,  this  is  a 
real-time,  state-space  model  of  the  system  to  be  controlled  that  predicts  its  behavior  and 
tracks  it  It  has  the  form: 
e{k)  =  y{k)-c-x(k ) 
x{k+\)  =  a-x{k)+b-u(k)+h-e(k ) 
u(k)  =  f  -x{k)  +  g-e(k) 

where  jc(&)  is  the  observer  state,  y{k)  is  the  set  of  sensor  signals,  e(k)  is  the  error  in 
predicting  the  output,  and  u(k)  is  the  actuator  command  vector.  The  matrices  a,  b ,  and 
c  are  identified  as  properties  of  the  system  to  be  controlled,  and  the  gains  /,  g,  and  h 
emerge  from  the  control  and  observer  optimization.  The  main  function  of  the  fast-rate 
controller  is  to  execute  these  equations  in  real  time.  These  equations  are  an 
implementation  of  the  optimal  LQG  control  law  for  the  identified  system.  Further 
implementation  details  are  provided  below  in  Section  0. 

In  parallel  with  the  Controller  is  the  Probe  Signal  Generator.  This  block  puts  colored 
noise  into  each  actuator  signal  to  use  as  a  persistent  excitation.  Coloring  is  accomplished 
by  means  of  a  comb  filter.  The  frequency  weighting  of  the  noise  is  calculated  in  the 
Probe  Signal  Gain  Adjuster,  based  on  the  accuracy  of  the  Controller’s  predictions  of  the 
sensor  signals  (see(8))  in  each  frequency  band.  Once  sufficient  accuracy  is  obtained  in 
any  band,  the  probe  signal  is  gradually  turned  off  in  that  band,  eventually  turning  off 
entirely  and  leaving  the  controller  to  operate  without  any  disturbance. 


The  most  important  block  in  the  operation  of  the  system  (other  than  the  Controller 
itself)  is  the  Plant  Identifier.  It  develops  a  model  of  both  the  effects  of  the  control 
actuators  on  the  sensor  signals  and  the  open-loop  disturbance  spectra  and  noises,  all  in  a 
state-space  model  form.  The  algorithm  is  Larimore’s  Canonical  Variate  Analysis 
method,  described  further  below. 

The  control  optimizer  is  an  ordinary  LQG  step  based  on  Riccati  equation  methods 
using  standard  Matlab  Riccati  solvers,  with  some  attention  paid  to  the  choices  of 
weighting  matrices  and  the  conditioning  of  inputs.  The  entire  set  of  software  is  driven 
from  Matlab.  The  organization  and  implementation  of  the  software,  and  its  relationship 
to  the  hardware,  is  described  below. 

3.3b.  Plant  Identifier  Theory  and  Implementation 

The  Plant  Identifier  (PI)  is  the  key  to  successful  operation  of  the  adaptive  controller. 

If  it  is  successful  in  characterizing  the  system  to  be  controlled,  then  the  controller  will 
certainly  work.  The  PI  attempts  to  identify,  using  only  time  histories  of  the  actuator 
commands  and  sensor  signals,  the  parameters  of  a  model  of  the  system.  It  uses  an 
implementation  of  Larimore’s  Canonical  Variate  method21'23  for  identifying  the 
parameters  of  a  state-space  model,  using  only  time  histories  of  the  exogenous  inputs  and 
outputs.  This  implementation  is  described  in  the  following  subsections. 

The  Problem  and  Top-Level  Solution  Description: 

The  identification  problem,  in  brief,  is  as  follows.  Given  a  set  of  time  histories  (possibly 
distinct  windows)  of  the  input  vectors  u{k) and  output  vectors  y{k) ,  a  model  is 
identified  of  the  form: 

x(k+\)  =  a-x(k)+b-u(k)+wx(k)  (11) 

y{k)  =  c-x{h)  +  d-u(k)+w2(k)  (12) 

where  x(k)  is  the  system  state,  and  (k)  and  w2  (k)  are  unknown,  zero-mean, 
inaccessible  white  noises  having  covariances  isjw1(£)w'1(&)rJ  =  ^, 

£|w2  (k)w2  (&)T|  =  V2 ,  and  Isjv^  (&)w2  (&)r  j  =  Vl2 .  The  model  parameters  identified 

are  thus  a ,  b ,  c,  Vx,  V2,  and  Vl2 .  The  matrix  d  may  be  included  or  not,  depending  on 
whether  a  feed  through  term  from  input  to  output  is  postulated  for  the  system  or  not.  The 
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method  works  by  finding  an  estimate  of  the  state  of  the  system,  and  then  finding  the 
parameters  of  the  model  by  least  squares  estimation,  assuming  the  state  is  as  estimated. 
The  order  of  the  system  is  chosen  by  evaluating  models  of  a  number  of  orders  and 
picking  the  “best”  one,  using  the  modified  Akaike  Information  Criterion  (AIC)24,25  as  the 
optimization  criterion. 

Three  essential  features  of  this  identification  algorithm  made  it  particularly  useful  for 
this  problem: 

1.  The  algorithm  works  in  the  presence  of  feedback.  This  is  shown,  for  example,  in 
Ref  22. 

2.  The  method  provides  an  effective  means  of  automatically  selecting  the  model 
order. 

3.  The  method  estimates  not  only  the  state-space  parameters,  but  the  statistics  of  the 
noise  as  well. 

Because  the  AIC  is  used  for  discriminating  among  various  possible  models,  we 
briefly  digress  to  describe  it’s  nature  qualitatively.  Our  self-contained  development  of 
the  theoiy  is  reserved  for  another  memo26.  The  AIC  is  an  information-theoretic  method 
for  evaluating  the  nearness  of  a  model  to  the  “true”  or  operating  one.  It  takes  into 
account  both  the  nearness  of  the  fit  of  the  model  to  the  data  and  the  number  of 
independent  parameters  that  were  adjusted  to  obtain  the  fit.  By  “model”  is  meant,  in  this 
context,  a  probability  density  function  for  the  observed  outputs  given  the  observed  inputs. 
The  AIC  works  for  any  model  in  which  the  history  of  outputs  is  a  deterministic 
(nonlinear)  function  of  the  inputs  and  the  model  parameters,  plus  an  additive  Gaussian 
noise,  provided  the  parameters  are  estimated  using  (an  approximation  to)  maximum 
likelihood  estimation  (MLE).  Recall  that  in  MLE,  one  selects  those  values  of  parameters 
that  maximize  the  predicted  probability  for  the  events  actually  observed.  Numerically, 
the  AIC  is  the  log  of  the  determinant  of  the  best  fit  error  covariance  plus  twice  the 
number  of  free  parameters  (an  additional  penalization  factor  is  applied  to  the  free 
parameters  for  small  samples).  The  absolute  level  of  the  AIC  is  not  significant,  only 
relative  changes,  because  we  are  comparing  a  number  of  models  in  order  to  choose  the 
best  one.  - 
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Returning  to  the  ID  procedure,  there  are  three  basic  steps: 

Decide  on  the  number  of  time  lags  Nlags  to  use  in  constructing  the  estimated  state. 
The  state  of  the  system,  in  keeping  with  linear  Markov  assumptions,  is  assumed  to  be  a 
linear  combination  of  past  inputs  and  outputs.  A  practical  decision  must  be  made  as  to 
how  many  such  past  outputs  will  be  considered  in  the  estimation  of  the  state.  This  is 
done  by  fitting  a  one-step-ahead  ARX  prediction  model  to  the  data  using  maximum 
likelihood  best  fits,  and  then  choosing  the  number  of  delays  that  produces  the  best-fit 
according  to  the  AIC.  The  user  must  specify  an  upper  bound  for  the  number  of  possible 
delays  to  consider. 

Estimate  the  state  history  for  the  system.  This  is  the  most  difficult  part  of  the  problem. 
According  to  Maikov  theory,  the  state  is  that  vector  that  is  just  as  useful  in  predicting  the 
future  outputs  as  is  the  entire  past  history  of  inputs  and  outputs.  The  CVA  procedure 
creates  a  future  vector  consisting  of  the  Nlags  next  output  values,  with  the  influence  of  all 
present  and future  input  values  removed  (present  value  influence  only  removed  if  there’s 
no  feedthrough).  To  remove  this  future  input  influence,  linear  regression  is  used.  The 
final  step  is  to  find  the  linear  combination  of  past  inputs  and  outputs  that  best  predicts  the 
future  vector,  for  each  possible  order  up  to  Nlags.  This  last  step  is  the  Canonical  Variates 
procedure,  and  is  accomplished  by  three  singular  value  decompositions. 

Estimate  model  parameters  for  each  model  order  and  choose  the  best  one.  The  model 
parameters  are  chosen  according  to  least  squares,  assuming  the  state  consists  of  the  first 
n  elements  of  the  state  history  vector  from  the  previous  step.  For  example,  to  estimate 
the  a  and  b  matrices,  those  values  of  matrix  are  chosen  that  minimize  the  error  in 
predicting  x(&  +  l)  given  x(k)  and  u(k),  as  can  be  seen  in(ll).  The  residual  error  is 

taken  as  an  estimate  of  the  white  disturbance  vector  (k) ,  which  is  then  used  to 
estimate  the  noise  covariances.  The  best  fit  error  for  the  model  is  not,  however,  just  the 
output  covariance  estimate  V2 ,  but  the  (larger)  Kalman  filter  prediction  covariance,  which 

depends  on  all  three  covariance  matrices*  as  well  as  the  other  parameters.  This 
covariance,  along  with  the  true  number  of  free  stochastic  model  parameters,  is  used  in  the 
evaluation  criterion  of  AIC.  All  possible  model  orders  are  evaluated,  and  the  best  is  then 
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selected  (i.e.,  that  minimizing  the  AIC),  and  its  model  parameters  recomputed  as  the  final 
estimates. 

These  basic  steps  are  elaborated  in  greater  detail  below,  once  we  elaborate  on  another 
method  of  Larimore’s27  for  saving  computations  in  least  squares  problems  in  which 
successively  more  regressor  variables  are  to  be  considered  in  the  fit. 

Order-Recursive  Linear  Least  Squares: 

Consider  the  traditional  problem  of  finding  the  matrix  M  that  minimizes  the  cost: 

J  =  tr[{Y-MX){Y-MX)T }  (13) 

In  this  problem,  the  columns  of  the  data  matrices  X  and  Y  might  represent  time  series. 

Suppose  we  have  a  matrix  'P  such  that  ¥ 7XX7y¥  =  Ip,  where  p  =  rank(X) .  Then 

XTx¥x¥TXXT  =  XT  as  can  be  verified  using  the  compact  SVD  for  X .  Expanding  the 
cost  we  get 

J  =  tr{YYT  -MXYT  -  YXTMT  +MXXTMT  } 

=  tr{YYT  -MXXt'¥'¥tXYt -YXt'¥'¥tXXtMt  +MXXtMt] 

{YY7 -YX7^7 XX7^7 XY7  1  (14) 

~tr\+YXT'¥'¥TXX7x¥'¥TXYT -MXXTxi,x¥TXYT  - YXT'¥'¥TXXTMT  +MXXTMT J 

=  tr[Y  (/  -  Xt'¥'¥tX)Yt)  +  tr^M  -  YXT'¥X¥T)XX7  (M  -  YXTx P'P7'  f  } 
yielding  a  global  minimum  of 

Jmia=tr{Y(l-XTxYxY7X)Y7}  (15) 

for  the  choice 

Mmm  =YX7'¥'¥t.  (16) 

Well,  suppose  we  have  another  set  of  regressor  variables  we  want  to  consider  in  addition 

to  X,  Z ,  and  suppose  we  have,  in  similar  fashion  to  before,  <hrZZr< I>  =  Iq,  where 
q  =  rank{Z } .  We  need  to  develop,  at  minimum  expense  of  computation,  a  new  matrix 


'F ,  such  that 
'X 


[XT  ZT]'Y  =  I-p, 


(17) 
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where  p  =  ra«&  j  ^  j .  The  new  least  squares  with  both  variables  as  regressors  is  then 
easy  to  compute,  as  before. 

The  method  is  as  follows.  Form  the  full  rectangular  SVD 
(<S>TZ)(XTx¥)  =  UCVT 

where  C  is  a  qxp  matrix,  and  the  singular  values  (diagonal  of  C )  satisfy  0  <  CiJ  <  1 . 

This  inequality  obtains  because  both  XTxVV  and  ZT<t>U  are  orthogonal  matrices,  with 
unit  vectors  in  the  columns,  so  that  C  represents  simply  the  mutual  dot  product  of  these 
unit  columns  of  the  two  matrices,  and  a  dot  product  of  unit  vectors  is  the  cosine  of  the 
angle  between  them.  Next,  suppose  the  number  of  diagonal  elements  of  C  of  value  less 
than  1  is  r .  Form 


Note  that 

T(I-CCt)Tt  =Ir. 

Finally,  set 

__[>F  -xVVCtTt 
O^n  ®UTT 

nzxp 

We  now  verify  that  this  matrix  satisfies(17): 


(19) 

(20) 

(21) 
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r  x  \xT  zr]'P  = 


VTx¥TXXT 

rTvfjT  \r\rT  ,  nrr 


-TCVTx¥T  TUt&t  j\_ZXT  ZZr J[0 
VTWrXZT 

T  7  VT  _jrr,l/TxUT  V7t  I  TTrT^T  77 T 


xxT  xzTl  w  -x vvcTr 


®UTT 


-TCVTxVT  XX  +  TU&  ZX  -TCVTx¥tXZt  +TUt®tZZt HO*,,  4>UTr 


yv  -'PFCrr 


VTx¥TXXT'¥V 


-VTx¥TXXTx¥VCTTT  +VTx¥tXZt®UTt 


-TCVTx¥t  XXTx¥V  +  TUt&tZXt  'F  V 


TCVTxVt  XX^VCT1  -TU1®1  ZXTx¥VCtT 


TCVTxVtXZt^UT1  +TU‘&ZZ,®UT 


T/fJ  77T, 


Ip  -CtTt+CtTt  1  - 

f TC'C'T  Tr  ^  z:  p 

-TC  +  TC  0  T(I-CCt)T 7 

\-tccttt  +nT  1  V  ' 


.0 


as  required. 

»  '  —  [X 

Intermediate  computations  are  made  easier  by  defining  X  =  'FrJT  and  X  =  X¥T 

Z 

“  '  A  ~ 

X1 

The  update  is  found  as  follows:  X'  -  VTX,  X  =  ,  „  «  .  . 

T(UZ-CX') 

Computations  of  this  type  are  found  in  both  the  first  and  second  major  steps  of  the 
algorithm. 

First  ID  Step:  ARX  Model  Estimation  to  Select  the  Number  of  Lags: 

An  ARX  model  is  of  the  form: 


mugs 

y(k)=  'EA,y(k-i)+'EB'u(k~i)+w(k)-  (22) 

i=l  i-0 

The  second  sum  starts  from  an  index  of  1  if  there  is  no  direct  feedthrough  postulated. 
The  procedure  is  to  compute  the  residual  error  from  a  linear  least  squares  fit  of  a  model 
having  a  number  of  lags  specified,  and  then  use  the  AIC  as  the  indicator  of  the  best  model 
order.  The  method  described  in  the  previous  section  is  used  to  gradually  add  more  past 
data  to  the  fit. 


Using  the  notation  of  the  previous  section  we  set: 
y  =  [y(0)  ...  y(tf-l)] 

x  \y{~m )  -  j'Ctf-m+i)! 


Z(m)  = 


u(-m )  •••  u(N-m  + 1) 
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where  Z(m)  means  the  Z  matrix  we’re  adding  to  investigate  m  lags.  Note  that  X  is 
defined  recursively: 

''»  =  Pri7~5l-  ^(0)  =  f“(0)  feed,Hm  (25) 

L  Z\m)  J  [  empty  else 

All  lag  orders  up  to  the  maximum  specified  by  the  user  are  examined,  and  the  one  giving 
the  least  AIC  is  used. 

Second  ID  Step:  State  History  Estimation: 

The  state  history  estimation  step  is  actually  comprised  of  three  principal  sub-steps: 
formation  of  the  past  history  vector,  formation  of  the  future  history  vector,  and  CVA 
transformation. 

The  first  sub-step  is  fairly  simple.  Considering  a  time  history  of  N  steps  beginning  with 
time  0,  the  past  history  matrix  is: 

y(-i)  - 

(— nlags )  •••  y(N-nlags) 

~u(-l)  ”  u{N-lj~ 

{-nlags)  •••  u(N -nlags) 

Column  ke  {0,1,...,jV-1}  represents  all  of  the  information  on  which  the  state  x  {k )  can 
depend. 

The  raw  future  matrix  is  simply: 

r  y{k) 

Fo=[iJ(0)  —  Hk)=  5  (27) 

y(k  +  nlags- 1) 

in  the  absence  of  any  external  input,  this  would  represent  all  of  the  future  output 
information  to  which  the  state  could  have  any  relevance.  In  the  normal  case,  when  input 
is  presence,  this  matrix  must  have  removed  from  it,  to  the  extent  possible  with  the 
information  available,  the  influence  of  future  inputs.  Consider  column  k  of  F0,  in  the 
case  when  no  feedthrough  is  assumed  in  the  plant.  The  influences  of  all  inputs  from 
u{k)  through  u{k  + nlags- 1)  must  be  removed.  Because  of  causality,  we  know 
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(without  feedthrough)  that  y{k)  can  only  be  affected  by  u(k- 1)  and  earlier  inputs. 
Ignoring  the  feedthrough  and  the  stochastic  components,  we  can  write: 

r~l 

y{k  +  r)  =  arx(k)  +  ^camb-u(k  +  r-m- 1).  (28) 

m=  0 

Define  the  sensor  reading  y(k + r)  without  the  previous  /  inputs  (noting  causality 
restrictions)  as: 

y(k  +  r)  i<  0 

/-I 

y(k+r)-'YJcamb-u{k+r-m-\)  0<i<r  (29) 

m-0 
r- 1 

y(k  +  r)-^camb-u(k  +  r-m-i)  i>r 

m-0 

The  future  vector  matrix  we  really  want  is: 

y(k,0,i) 

F  =  [y(0,nlags-l)  •••  y(N-\,nlags-\)~],  y(k,i)=  y(k,r,i)  (30) 

y(k,nlags-l,i) 

The  basic  idea  is,  we  start  with  the  raw  future  vector,  which  is  j>(£,0) ,  and  iterate  it 
forward  to  y(k,nlags-\) ,  which  is  the  desired  future  vector.  The  key  iteration  is: 

y(k,r,i+i)  =  y(k,r,i)-ca'b-u(k  +  r— i— 1),  i<r,  (31) 

which  only  need  by  performed  on  the  bottom  part  of  the  vector  each  time  (one  fewer  each 
iteration).  But  whence  ca'b  ?  Plug  in  to  get  the  explicit  expression: 

y(k, / + 1,  /)  =  aMx(k)+ calb  -u(k).  (32) 

We  don’t  have  the  state — that’s  the  whole  idea,  after  all — but  suppose  we  have  some 
vector  t|  (&)  containing  the  state: 

x(*)  =  £il(jfc)  (33) 

Then 

y(*,/+U) =[cJb  «“£]  .  (34) 

We  could  get  this  matrix  by  regression,  and  then  just  pull  out  the  first  nu  columns  to  get 
ca'b : 
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c-Hm-'Ho ...  y(*-u+u)][2j : 

But  notice  the  product  of  the  last  two  matrices  is  independent  of  i ,  so  we  can  hold  onto 


uo)  -  u(w-i)Trc 
L'n(o)  -  ti(at-i)J  L o 

For  the  t|(Ar)  vector  we  just  use  the  column  space  basis  vector  obtained  from  an  SVD  of 
the  past  matrix  from  (26). 

The  summary  of  the  method  for  obtaining  the  future  vector  free  of  future  u  is  thus: 

1.  Set  up  the  initial,  raw,  future  matrix  as  in  (27).  Compute  T  as  in  (36).  Set  /= 0 

2.  Compute  ca'b  via  (35).  Subtract  out  the  influence  of  u(k  +  r-i- 1)  from  row 
r+1  of  F ,  for  r  =  i+l:nlags-l.  Quit  if  i=nlags- 2,  otherwise  increment  i 

f  ‘ 

and  repeat. 

The  final  state  estimation  step  is  to  perform  a  CVA  relating  the  past  to  the  future  matrix. 
First  we  find,  via  compact  (square  inner  matrix)  SVD: 

P  =  Up (37) 
and  similarly  for  the  future  vector: 

F  =  U/LfV/.  (38) 

Note  that,  for  the  purposes  of  the  information  present  in  the  time  history  vectors,  the 
column  space  bases  Vp  and  Vf  are  excellent,  perfectly  conditioned  (e.g.  VpV  =  I p ) 

substitutes  for  the  originals.  The  relationship  between  them  is  revealed  when  we 
compute  the  full  SVD  of  the  dot-product  matrix: 

VfTV,=UJLjrJ  (39) 

The  state  is  simply 

X  =  Vj-VpT.  (40) 

This  produces  a  set  of  states  that  are  independent  of  one  another,  related  one-to-one  with 
future  information  vectors,  and  ordered  in  descending  order  of  importance. 

More  can  be  found  on  the  CVA  and  its  relevance  to  MLE  in  Ref.  26. 
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Third  ID  Step:  Estimate  Model  Parameters  and  Pick  the  Best  One 
The  last  step  in  the  ID  procedure  is  fairly  simple.  The  parameters  of  the  model  are 
computed  using  least  squares,  for  each  possible  model  order  up  to  the  size  of  the  system 
state  calculated  in  the  previous  step.  When  there  is  no  feedthrough  postulated,  two  least 
squares  problems  emerge.  First,  the  state  update  equation  is  used.  The  state  update 
equation  can  be  assembled  as: 

0(1)  -  *(*)]  =  [»  a]“{o)  I  "[tfllj  +[”'.(°)  -  ».(*-<)]  (41) 

where  the  state  history  and  the  control  input  history  are  given,  the  a  and  b  matrices  are 
to  be  estimated,  and  the  noise  history  w,  (A)  is  taken  as  the  least-squares  fit  error.  If  a 

model  of  order  n  is  being  estimated,  then  only  the  first  n  state  values,  as  calculated  in 
the  previous  major  step,  are  included  in  the  fit.  The  second  fit  problem  uses  the  measured 
sensor  readings: 

[y(0)  •••  y(W-l)]  =  c[x(0)  —  x(JV-l)]  +  [w2(0)  •••  w2(W-l)]  (42) 

where  now  the  noise  history  w2  (&)  is  taken  as  the  fit  error. 

The  covariances  are  estimated  via  the  residuals  using  simple  averaging,  e.g., 

^-^■[»,(°)  -  '•’,(*-!)][>. (°)  -  wI(JV-l)]r.  (43) 

I 

The  case  of  feedthrough  is  somewhat  simpler,  because  only  one,  combined,  least  squares 
problem  arises: 

x(l)  •••  x(N)  _  b  a  m(0)  •••  tt(N-l)  w,(0)  •••  w^W-l) 

y(0)  •••  ^(W-l)  d  c  x(0)  •••  x(iV-l)  +  w2(0)  •••  w2(N-1)  ' 

These  least  squares  problems  are  managed  as  in  the  section  on  re  to  be  considered  in  the 
fit. 

Order-Recursive  Linear  Least  Squares,  so  that  the  inclusion  of  the  next  higher  order  state 
in  the  problem  has  minimal  computational  impact. 

An  important  subtlety  is  the  assessment  of  the  AIC  for  each  model  order.  This  point  is 
left  (intentionally?)  vague  in  Larimore’s  papers.  It  seems  clear  that  the  simple  average  fit 
error  V2  should  not  be  used,  because  the  errors  from  different  steps  are  not  independent, 
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so  that  the  standard  independent  trials  arguments  of  MLE  do  not  apply.  Instead,  the 
innovations  process  of  the  associated  Kalman  filter,  which  in  the  steady  state  produces 
independent,  identically  distributed  errors,  must  be  used.  This  is  computed  via  the 
Algebraic  Riccati  Equation  (ARE)  of  Kalman  filtering.  First,  find  q  to  solve: 

q  =  a-q-aT  +Vv-{a-q-cT  +Vn)(c-q-cT  ^V2Y(a-q-cT  +  (45) 

Then,  the  innovations  covariance  is: 

V  =  c-q  cT  +V2  (46) 

In  finding  the  solution  to  the  ARE,  it  is  necessary  to  ensure  that  the  states  are  detectable 
(stable  or  observable) — which  may  not  be  true  numerically  if  the  order  is  high. 
Preconditioning  is  used  to  ensure  both  stability  and  observeability  in  our  implementation. 
Another  method  for  obtaining  V  would  be  to  run  the  Kalman  filter  to  get  the  innovations 
process,  and  then  use  the  innovations  to  obtain  essentially  the  same  result  as  (44).  To  get 
the  Kalman  filter  would  require  the  same  ARE  solution  procedure. 

The  AIC  is  then  computed  for  each  model  order  as: 


J4/C  =  iV(ny(l  +  log(27t)))  +  log(det(F))  +  2/  m  (47) 

where  N  is  the  number  of  output  time  steps  examined,  and  the  number  of  parameters  m 
is 


m  =  < 


—ny(ny  +  \)  +  (2ny  +  nu)n  d  =  0 

2 

— ny{ny  +  \)  +  (2ny  +  nu)n  +  ny-nu  else 


(48) 


The  number  of  states  is  «,  the  number  of  measurements  ny ,  the  number  of  actuation 
inputs  nu .  Notice  that  this  number  of  parameters  is  far  fewer  than  the  number  of 
elements  in  the  matrices  of  the  identified  model,  because  there  are  a  number  of  free 
parameters  still  remaining,  which  are  arbitrary.  The  number  of  fixed  parameters  was 
derived  in  Ref.  31.  and  the  stochastic  part  of  this  derivation  can  be  also  be  done  by 
placing  the  plant  in  an  equivalent  Kalman-filter  form: 
x(£+l)  =  a-;c(A)+6-w(&)+/j-w(&) 


y(£)  =  c-x(£)+w(£),  £’|w(A)w(A:)rJ  =  F , 


(49) 

(50) 
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where  the  h  and  V  matrices  are  from  the  Kalman  filter  computation  as  before.  It  can  be 
shown  that  this  rendering  of  the  plant  model  (not  the  controller!)  has  exactly  the  same 
statistical  behavior  (output  correlation  function)  as  the  original  plant  model.  The  total 
number  of  parameters  is  then  that  of  the  symmetric  V  {ny{ny  +  \)!2\  of  h  (n-ny),  plus 


that  of  the  a-b-c  realization  ( n  •  [ny +nu ) ),  keeping  in  mind  the  remaining  free  state  basis 

change  allowed,  which  eliminates  n2  parameters  from  the  bare  matrix  count. 

The  factor  /  provides  an  extra  penalty  on  the  number  of  parameters  when  the  sample 
size  is  small.  It  is  calculated  as: 

N 


/  = 


N -—-^-ny{ny  +  \) 
ny  2 


(51) 


After  selecting  the  best  model  order,  the  parameters  (having  previously  been  discarded) 
are  then  recalculated,  completing  the  identification  process. 


3.4.  Control  Algorithm  Theory  and  Implementation 

The  control  algorithm  uses  at  its  core  the  Linear  Quadratic  Gaussian  theory 
available  since  the  1960s30.  This  theory  solves  the  optimum  least-squares,  steady-state, 
infinite  horizon  control  problem  for  a  linear  time-invariant  system.  To  facilitate 
discussion,  we  review  the  theory  results  very  briefly  (to  the  extent  of  restating  the 
complete  design  equations,  and  adapting  them  to  our  implementation),  and  then  proceed 
to  a  discussion  of  some  of  the  practical  aspects  of  the  implementation  (for  example,  how 
do  we  update  the  controller  on  the  fly?). 

LQG  Theory  Summary: 

Despite  the  ubiquity  of  LQG  theory,  it  is  remarkably  difficult  to  find  a  really 
complete  set  of  design  equations  written  down  anywhere.  Reference  3 1  applies  to  the 
more  difficult  problem  of  designing  an  optimum  order-constrained  controller  for  a 
higher-order  system.  However,  the  equations  are  easily  reduced  to  the  LQG  situation  as  a 
special  case:  we  present  the  results  of  this  reduction  here.  The  problem  is  to  minimize, 
by  means  of  dynamic  output  feedback  compensation  applied  to  the  system  we  described 
earlier  in  (1 1)  and  (12)  (no  feedthrough  term),  the  cost  functional: 
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The  penalty  is  used  to  control  actuator  effort,  and  the  coupling  matrix  Rl2  is  zero  in 


most  problems.  The  controller  is  constrained  to  be  of  the  form: 

xc{k  +  \)  =  ac-xc{k)+bc-y{k)  (53) 

u(k)  =  cc-xc{k)  +  dc-y(k).  (54) 

The  design  procedure,  given  these  data,  is  to  solve  the  following  two  ARE’s  for  positive 
semi-definite  p  and  q : 

q  =  a ■q-aT +Vl-(a-q-cT +Vl2)(c-q-cT +  V2)  '  (a-q-cT +Vl2f  (55) 

p  =  aT  •  p-a  +  R^  ~{aT  •  p -b  +  /^2)(&r  -  p  b  +  1  (ar  •  p •  b  +  /?12)  (56) 

The  optimal  matrix  values  are  then: 

dc  =  -(bT -p-b  +  Rtf1  (bT -p-a-q-cT  +RnT -q-cT  +bT -p-V^cq-c7  +V2y  (57) 
bc=(a-q-cT  +Vu)(cq-cT  +V2)1  +b-dc  (58) 

ce=-(bT  ■pb  +  R1)\aT -pb  +  Ruf  -dc-c  (59) 

ac=a+b-dc-c-bc-c+bcc  (60) 


The  dynamic  compensation  form  of  (53)  and  (54)  does  not  bring  out  the  observer  that  is 
embedded  in  the  controller,  which  we  need  to  generate  and  monitor  the  modeling  error  in 
real  time.  The  observer  from  we  have  implemented  is  developed  as  follows.  Define: 


h  =  (aq-cT +Vn)(c-q-cT  +  V2)  (61) 

/  =  (62) 

g  =  dc  (63) 

Then  we  have: 

bc  =  h  +  b-g  (64) 
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Cc=f-g-c  (65) 

ac=a  +  b-g-c-(h  +  bg)c  +  b-(f-g-c)  =  a-bg-c-h-c  +  b-f  (66) 

So  that 

u{k)  =  {f-g-c)xc(k)  +  g-y(k)  =  f-xc(k)  +  g(y(k)-c-xc(k))  (67) 

Letting 

e{k)  =  y(k)-c-xc(k)  (68) 

be  the  output  prediction  error  (i.e.,  the  innovation),  the  control  becomes 

u(k)  =  f-xc(k)  +  g-e(k).  (69) 

Turning  now  to  the  state  update, 

xc(k  +  l)  =  (a-b-g-c-h-c  +  b-f)xc(k)  +  (h  +  b-g)y(k) 

-a-xc(k)  +  b[f-xc(k)  +  g(y(k)-c-xc(k)))  +  h{y(k)-c-xc(k))  (70) 

=  a-xc(k)+b-u(k)  +  h-e(k) 


so  that  the  state  update  is  just  the  same  as  the  plant  with  Kalman  observer  innovation 
corrections.  The  control  is  plant  state  estimate  feedback  with  innovation  direct 
feedthrough.  This  rendering  of  the  control  system  observer  is  the  method  of 
implementation  in  the  code. 

Performance  Monitoring  and  Control  Reversion: 

It  is  not  unlikely  that  at  some  point  during  the  convergence  process,  the  inaccuracies 
of  the  identified  plant  relative  to  the  true  system  will  be  such  as  to  cause  the  optimal 
controller  for  the  estimated  plant  to  be  destabilizing  when  connected  with  the  true  system. 
To  guard  against  this  eventuality,  each  new  controller  is  tested  on  a  trial  basis  for  a  pre¬ 
specified  amount  of  time  before  being  accepted  as  the  new  baseline.  If,  during  the  trial 
time,  the  new  controller  is  deemed  unacceptable,  the  real-time  system  reverts  back  to  the 
previous  baseline.  The  specific  sequence  is  as  follows,  once  a  new  controller  is  available 
for  upload: 

1.  Compute  matrix  transformations  to  transform  the  state  of  the  baseline  controller 
to  that  of  the  new  controller,  and  back  again  (in  case  reversion  is  necessary). 

2.  Upload  the  new  controller  and  associated  conversion  matrices. 
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3.  Calculate  the  new  controller  state,  start  the  trial  timer,  reset  the  performance 
monitors,  and  begin  executing  the  new  controller. 

4.  If  the  performance  monitors  detect  a  fault,  map  the  new  controller  state  back  to 
the  baseline  controller  state,  and  revert  to  the  baseline  (raising  appropriate  fault 
flags). 

5.  If  the  trial  time  runs  out,  without  any  performance  faults,  accept  the  new 
controller  as  the  baseline  and  set  appropriate  flags. 

There  are  three  performance  monitors  that  can  raise  an  objection  to  the  new  controller 
actuator  effort,  output  prediction  error,  and  sensor  response.  These  monitors  function  by 
simply  filtering  the  square  of  each  actuator,  sensor,  or  prediction  error  variable,  as 
follows: 

?.(0=f.(0+8,(®(02-«.(0)  (71) 


where  qe  (/)  is  the  running  average  of  the  i,h  prediction  error.  If  any  of  these  running 

averages  exceeds  the  threshold  set  for  it  during  the  trial  time,  the  controller  is  deemed 
unacceptable. 

It  remains  for  us  to  clarify  how  the  state  of  the  baseline  observer  is  converted  to  an 
“equivalent”  state  in  the  new,  trial  one.  We  choose  that  state  in  the  new  observer  that, 
given  no  further  control  inputs,  produces  the  mean-square  closest  history  of  sensor 
signals  to  that  of  the  baseline.  That  is,  we  choose  x2  to  minimize 


qo  2 

J  -UK*,  ~c2a2kx2 


=  XlT  k  ~2XlT  YualkciC2a2  x2  +  X2  X 


Consider  the  following  iteration: 


rc,(*+i)  Q^+OVlY  aw  aiWYV  oYfV  v/1  (73) 

Qn  (k  +  1)  Q2(k  + 1),  (o  a2j  ,Q\2^)  Ql  (k)  ^  0  ai  i  \C2CI  C2C2  2 

with 
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0i  (°)  Qn  (°)1  _(  cxc ,r  -c,c2 

Qn  (0)  02  (0);  l-C2Cl  C2C2  ) 


It’s  not  hard  to  verify: 

Qi  (*)  =  Qu  (k)  =  Yu^Cx^al,  Q2  (k)  = 


The  solution  to  the  Algebraic  Lyapunov  Equation  produces  the  limiting  sum  as  k  -» °° : 

Ql  Q*)J°1  °¥fl 

^12T  02  J  1°  aJlfi/  02  A  0  a/J  lc2c,r  c2c2r  / 


J  —  X [  Q\X\  2xi  0X2^2  *2  02*2  • 

The  optimum  choice  is  then  found  by  rearrangement 

J  =  xi  Q\XX  —  xi  01202  02*2  —  *2  0202  012  *1  +  *2  02*2 

=  *ir0i*l— *1  01202  012  *1 

"*"*1  01202  0202  012  *1  _  *1  01202  02*2  —  *2  0202  012  *1  *2  02*2 
=  Xl  (01  —  01202  012  )^1— (*2  —  02  012  *l)  02  (*2  —  02  012  *1  ) 


i-e., 

*2=02+0!2r*l-  (79) 

So  the  conversion  matrix  is  02+0I2r »  and  in  like  fashion,  the  conversion  back  is  QxQn  • 

Probe  Signal  Theory  and  Implementation: 

The  probe  signal  is  used  before  proper  convergence  has  been  achieved  to  ensure  that 
sufficiently  accurate  system  identification  has  occurred  for  high  performance  in  each  sub¬ 
band  examined. 

The  probe  signal  vector  is  given  by: 

X(*)=£y,Re{M*)}  (80) 

F=Pl 
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where  %(&)  is  the  probe  signal  vector,  Bp(k)  is  the  complex  comb  filter  state  for 
frequency  bucket  p,  yp  is  the  real  scalar  gain  to  use  for  frequency  bucket  p ,  and  pl 
and  p2  are  the  indices  of  the  first  and  last  frequency  buckets  to  use. 

The  comb  filter  state  is  a  running  FFT  over  the  last  L  values  of  a  zero-mean  univariate 
white  noise  input  vector  w( k) : 


M*)= 


2  £4  n  .x  Mn  L 

-YJw{k-i)e  L  1  <P<- 

L  (=o  ^ 

I  L- 1  ,2*  , 

~y'<w{k~i)e  L  else 
L  i= o 


The  factor  of  two  is  used  for  all  but  the  DC  and  LI 2  frequency  buckets  because  only 
positive  frequency  buckets  are  included.  The  updates  to  6p  (£)  are  computed  recursively 
via: 

e'T"  I <p<^ 

M*H  »  ,  •  .  <82> 

e1  L?  •§p{k-\)+—[w(k)-w(k.-L))  else 

A  circular  buffer  is  used  to  maintain  a  record  of  the  last  L  values  of  the  input  noise. 

The  gains  yp  are  adjusted  as  follows: 

1.  Obtain  a  set  of  time  data  of  the  output  prediction  error.  Organize  the  data  into 
windows  of  length  L .  Compute  the  FFT  of  the  data. 

2.  In  each  frequency  bucket,  if  any  component  of  the  prediction  error  in  the  bucket  is 
larger  than  the  reference  level  for  that  component,  double  the  probe  gain  for  that 
bucket.  Otherwise,  cut  the  probe  gain  for  that  frequency  bucket  in  half. 

3.  Limit  the  allowable  probe  gain  level  in  each  frequency  bucket  to  a  predetermined 
amount. 

4.  Update  the  gains  in  the  probe  output  computation. 

This  completes  the  probe  computation  and  adjustment  theory. 
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Software  and  Hardware  Implementation: 

The  control  algorithm  is  implemented  using  a  combination  of  DSP-based  code  for 
fast  execution  time,  and  Matlab  routines  for  high  flexibility  and  ease  of  programming. 

The  block  diagram  of  the  hardware  implementation  is  shown  below  in  Figure  30.  A 
Gateway  PC  Workstation  hosts  an  Innovative  Integration  M67  DSP  board  in  one  of  its 
PCI  slots.  The  DSP  board  has  a  Texas  Instruments  ‘C67  digital  signal  processor,  and 
includes  a  Servo  16  analog  I/O  module  (also  supplies  by  Innovative  Integration).  The 
analog  signals  (both  input  and  output)  are  filtered  using  the  Krohn-Hite  anti-aliasing 
filtering  electronics  (comer  frequency  fixed  at  5000  Hz).  The  filtered  actuator  command 
voltage  is  then  routed  to  the  piezo  amplifier  (RCV  1000/7  made  by  Piezo-Mechanik  of 
Germany).  The  (0  to  1000  volt)  drive  voltage  from  this  amplifier  goes  to  the  piezo¬ 
actuator  (two  different  model  numbers  “PSt  1000/10/150  VS25  -HP”,  and  “PSt 
1000/10/150  VS25  -HS/HAT”  were  purchased,  also  from  Piezo-Mechanik,  to  compare 
their  stroke  at  high  frequency). 


Figure  30.  Block  Diagram  of  System  Hardware 

As  mentioned  previously,  the  algorithm  was  partitioned  to  executed  on  both  the 
DSP  and  the  PC,  as  shown  in  Figure  31  below.  Routines  requiring  fast  execution  time, 
such  as  data  acquisition  and  controller  implementation,  were  implemented  on  the  DSP. 
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Routines  that  could  be  run  in  the  background,  such  as  identification  and  control  design, 
were  run  in  Matlab  on  the  PC  (this  turned  out  to  be  huge  time  saver  because  it  eliminated 
the  conversion  of  Matlab  m-files  to  c-code).  Of  course,  other  routines  were  necessary  to 
transfer  commands  and  data  between  the  DSP  and  PC. 


Figure  31 .  Top  Level  Software  Architecture 


The  routines  of  highest  interest  working  on  the  PC  Host  in  the  Matlab  environment  are 
described  below  in  Table  3.  These  files  (with  the  exception  of  runTargetm)  are  in 
“c:\projects\common  mfiles”  The  top  level  routines  are  usually  called  directly  by  a 
program  such  as  “test3.m”.  In  the  table,  the  indented  routines  are  called  by  the  one  above 
(“adaptcontroLm”). 
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Table  3.  Main  Matlab  m-Files  for  Control  on  the  Host 


Matlab  m-File 

Purpose 

runTarget 

initializes  DSP  hardware,  starts  main  shell  and  interrupts 

adaptparams 

sets  up  default  parameters  for  the  adaptive  controller 

adaptcontrol 

main  routine  for  adaptive  control  on  host 

sendcmd 

sends  a  command  to  the  DSP  application 

getctlstatus 

gets  status  of  the  currently  operating  controller  from  the  DSP 

probeadjust 

finds  new  adjusted  gains  of  the  probe  frequency  buckets 

sendprobegains 

sends  the  host  the  new  probe  gains  for  immediate  implementation 

gethistory 

gets  another  time  window  of  data  from  the  DSP  of  actuator 

commands,  sensor  signals,  and  output  prediction  errors 

cvaid 

identifies  a  new  plant  model 

dlqg 

computes  the  parameters  of  an  LQG  controller 

stateconvert 

computes  conversion  of  states  between  two  controllers 

tiyobserver 

sends  a  new  controller  to  the  DSP  for  immediate  implementation, 

waits  to  see  whether  it  worked  out 

A  description  of  the  low-level  Matlab  m-files  and  related  files  for  controlling  and 
communicating  with  the  DSP  is  given  in  Table  4  below. 
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Table  4.  Matlab  m-Files  for  Use  in  DSP  Communication/Configuration 


file 

location  on  PC 

description 

m-files 

runTarget.m 

c:\MatlabRl  2\toolbox\PSI 

main  routine  for  running  an  application  on  the  DSP 

board 

genRTS  tructure.  m 

builds  the  real-time  (RT)  structure  used  to  download 

an  application  to  die  DSP  board 

initTargetm 

routine  that  initializes  communications  with  and 

downloads  the  application  to  the  DSP  board 

RTConfiglnitialize.m 

routine  that  configures  the  ‘real-time5  info:  sample 

rate,  active  A/D  and  D/A  channels,  etc. 

RTConfigDisplay.m 

routine  that  displays  ‘real-time5  configuration  to 

Matlab  workspace  (as  read  back  from  DSP  board) 

TransConfigInitialize.m 

routine  that  configures  the  ‘transfer5  info:  shared 

memory  location,  size 

TransConfigDisplay.m 

routine  that  displays  ‘transfer5  configuration  to 

Matlab  workspace  (as  read  back  from  DSP  board) 

DebugDisplay.m 

reads  and  displays  status  memory  from  ‘shared5 

memory  area 

support  files 

IlboardXControl  1  .ocx 

c:\M6x\DspComp\ActiveX 

ActiveX  control  for  communicating  with  M67  DSP 

board 

Two  main  DSP  C-code  applications  were  developed  on  this  project:  One  for  simply 
acquiring  and  storing  data  from  the  DSP  on  the  host,  and  the  other  for  implementing  the 
adaptive  control  algorithm.  A  description  of  the  DSP  C-code  software  files  is  given 
below  in  Table  5.  All  of  the  DSP  code  was  developed  using  Texas  Instruments  Code 
Composer  (Version  2.0). 
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Table  5.  C-Code  Software  Files  for  Use  on  the  DSP 


file 

location  on  PC 

description 

Data  Acquisition  Application 

dataAcq.pjt 

c:\projects\dataAcq 

Code  Composer  project  file 

dataAcq.c 

main  routine  for  implementing  DSP  data 

acquisition  routine 

Adaptive  Control  Application 

adapt  l.pjt 

c:\projects\adapt- 

M6x  code 

Code  Composer  project  file 

adaptMain.c 

main  routine  for  implementing  DSP  adaptive 

control  algorithm 

adapt_TI.c 

transferring  and  implementing  adaptive  routines 

observer_TI.c 

transferring  and  implementing  observer/controller 

probe_TI.c 

transferring  and  implementing  probe  signal 

comb_TI.c 

transferring  and  implementing  comb  filter 

realbuf_TI.c 

routines  for  capturing  stream  of  data 

plant.c 

example  plant  model 

I 

Support  Files  and  Libraries 

matm.c,.lib 

c:\C  Libraries 

misc.  matrix  and  vector  operations 

hosttalk.c,  .lib 

c:\C  Libraries 

higher  level  routines  for  processing  commands, 

and  transferring  data 

MIMOState_Space 

.c,  .lib 

c:\C  Libraries 

routines  for  implementing  a  MIMO  state  space 

controller 

realtime_M6x.c, 

.lib 

c:\C  Libraries 

routines  that  set-up  and  support  data  acquisition 

on  Servo  1 6  and  A4D4  I/O  modules 

transfer.c,  .lib 

c:\C  Libraries 

routines  for  transfer  data  to/from  host  PC 

generic.cmd 

c:\M6x 

memory  map  file 

3.5.  Control  Demonstration  Description  and  Results 

A  colored  random  noise  cancellation  demonstration  was  performed  in  PSI’s 
laboratory  space  to  demonstrate  the  operation  of  the  algorithm.  The  basic  control  setup 
was  as  shown  in  Figure  32  below. 

Open  Lab  Environment 


< -  70  inches  — - ► 


Figure  32.  Layout  of  the  Cancellation  Experiment 

A  speaker  was  used  as  a  disturbance  and  cancellation  source  for  a  random  noise 
cancellation  experiment.  Sound  was  measured  by  means  of  a  nearby  microphone  (used 
for  feedback  cancellation),  and  performance  was  monitored  approximately  70  inches 
away  by  a  second  microphone.  The  environment  was  an  open  laboratory  with  a  nearby 
wall,  door,  workbenches,  etc..  The  floor  was  carpeted. 

The  noise  source  was  generated  by  a  Siglab  2-channel  analyzer  connected  to  a 
lap-top  computer  running  Matlab.  Both  microphones  were  monitored  at  the  Siglab.  The 
noise  was  random  between  0-200  Hz  with  an  internal  high-order  low-pass  filter.  Both  the 
feedback  microphone  signal  and  the  noise  were  fed  to  the  control  DSP,  which  used  only 


61 


the  microphone  signal  for  control  adaptation.  The  noise  was  simply  added  to  the  control 
signal  drive. 

A  variety  of  sample  rates  were  used  for  experimentation  in  the  DSP,  with  the 
most  successful  experiments  run  at  the  highest  rate  (14  kHz).  The  results  of  such  an 
experiment,  with  an  observer  of  40  states  maximum,  are  shown  below  in  Figure  33. 
Approximately  10  dB  of  attenuation  was  achieved  over  the  significant  100-200  Hz  band. 
While  experimentation  showed  that  multi-path  and  3-D  transmission  effects  were 
significant,  performance  was  almost  certainly  limited  chiefly  by  time-delays  through  the 
drive  amplifier,  processor,  and  speaker. 

Convergence  was  quite  rapid  (a  few  time-windows  of  4096  points  of  data  at  14 
kHz)  for  the  observer  model.  Controller  design  was  the  limiting  factor,  because  a  slow 
loop  was  used  for  gradually  increasing  the  controller  authority  (this  can  be  modified,  but 
was  not,  due  to  time  constraints). 


io2  id3 

Frequency  (Hz) 


Figure  33.  Far-Field  Cancellation  Results  Using  Adaptive  Algorithm 
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It  was  necessary  artificially  to  increase  the  measurement  noise  parameter 
identified  by  the  system  prior  to  use  in  adaptive  design  to  ensure  that  the  resulting 
controller  was  internally  stable.  As  is  common  with  LQG,  too  high  a  control  effort  under 
low-noise  conditions  tends  to  cause  the  controller  to  be  very  aggressive,  introducing 
unstable  poles  that,  while  ostensibly  closed-loop  stable  when  connected  to  the  system, 
make  the  controller  extremely  sensitive  to  modeling  inaccuracies.  Overall,  the 
performance  seemed  close  to  the  physical  limits  of  the  laboratory  setup,  which  seemed 
dominated  by  the  lags  and  time  delays,  as  mentioned  previously. 


63 


4.  PULSE  MICROJET  ACTUATOR 

High  frequency  forcing  was  one  of  the  numerous  concepts  that  was  proposed  for 
reducing  high-speed  jet  noise.  The  goal  of  the  technique  is  to  affect  the  far-field  jet  noise 
radiation  by  affecting  the  airflow  before  it  exits  the  nozzle.  In  this  case,  high-pressure 
microjets  were  employed  to  affect  the  airflow  upstream  of  the  nozzle  exit.  The  microjets 
were  operated  at  high  frequencies  (2  kHz)  in  order  to  provide  an  unsteady  disturbance  to 
the  airflow. 

Prior  to  conducting  any  experiments,  it  was  necessary  to  design  and  fabricate  a 
nozzle  that  would  allow  for  the  microjet  injection  into  the  airflow  upstream  of  the  nozzle 
exit.  A  picture  of  the  nozzle  that  was  used  for  the  experiment  is  shown  in  Figure  34.  The 
converging  nozzle  has  an  exit  diameter  of  50.8  mm.  The  nozzle  had  four  ports  for 
microjet  injection  and  it  also  had  four  pressure  ports  before  and  after  the  microjets.  This 
allowed  us  to  measure  the  static  pressure  before  and  after  the  microjet  injection  so  that 
we  could  determine  changes  brought  about  by  the  microjet  injection  on  the  airflow.  The 
nozzle  was  constructed  of  Nistele  230,  a  nickel-based  material  that  allowed  us  to  operate 
the  jet  at  high  stagnation  temperatures  without  any  active  cooling. 


Figure  34.  Nozzle  for  High  Frequency  Forcing 
Experiments. 

It  was  also  necessary  to  procure  a  valve  that  could  handle  high  pressures  and  operate  at 
frequencies  up  to  2  kHz  without  any  active  cooling.  Such  a  valve  was  not  found  to  be 
commercially  available.  Therefore,  Planning  Systems  Incorporated  (PSI)  was  contracted 
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to  develop  and  provide  a  reasonably  priced  high-speed  valve  that  could  handle  the  high 
pressures.  The  intent  was  to  then  purchase  4  of  these  valves  for  use  on  the  four  ports  on 
the  nozzle  discussed  above.  A  picture  of  the  valve  that  was  provided  to  us  at  the  end  of 


the  contract  is  found  in  Figure  35.  This  valve  was  expensive,  difficult  to  operate  and  it 
required  active  cooling.  As  a  result,  a  new  high-speed  valve  was  developed  to  meet  the 
requirements  at  FMRL  (Fluid  Mechanics  Research  Laboratory).  A  schematic  of  the  valve 


Figure  36.  Simple  schematic  of  the  high-speed  valve  developed 
at  FMRL. 


Is  shown  in  Figure  36.  This  valve  easily  met  all  of  the  stated  requirements.  The  only 
drawback  was  that  it  required  a  high  torque  motor  in  order  to  operate  at  high  pressures. 
Preliminary  results  obtained  for  the  valve  can  be  found  in  Figure  36  and  37.  Figure  36  is 
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the  pressure-time  plot  for  the  valve  output  when  500  psi  is  applied  to  the  valve.  The 
pressure  was  measured  at  a  distance  of  5  diameters  downstream  of  the  valve  exit.  The 
mean  pressure  at  this  point  is  42.5  psi  with  fluctuations  of  17.5  psi.  Figure  37  is  the 
corresponding  narrowband  frequency  spectrum.  As  can  be  seen,  the  valve  is  operating 
well  above  the  required  2  kHz.  The  other  noise  peaks  in  the  plot,  except  for  the  harmonic 
at  4.5  kHz,  are  believed  to  be  a  result  of  the  high  torque  motor.  Increasing  or  decreasing 
the  pressure  applied  to  it  could  change  the  operational  speed  of  the  valve.  An  increase  in 
pressure  caused  the  valve  to  slow  down  and  decrease  in  pressure  caused  it  to  speed  up. 


Figure  36.  Pressure-time  plot  for  the  Figure  37.  Narrowband  frequency 

high-speed  valve  operated  at  500  psi  spectra  corresponding  to  Figure  36. 


Once  the  nozzle  and  the  valve  were  completed  it  was  possible  to  conduct  the 
experiment  investigating  the  effects  of  high  frequency  forcing.  A  picture  of  the  setup  for 
this  experiment  can  be  found  in  Figure  38.  The  high-speed  valve  was  used  to  supply  all 
four  of  the  microjets  on  the  nozzle.  The  microjets  were  operated  at  2.2  kHz  and  a 
stagnation  pressure  of  500  psi.  The  main  jet  was  operated  at  a  stagnation  pressure  of  44 
psi  and  a  stagnation  temperature  of  1033  K.  The  far-field  noise  radiation  was  acquired 
through  the  use  of  8  Bruel  &  Kjaer  4939  microphones.  These  microphones  were  set  up  at 
a  distance  of  65  nozzle  diameters  in  a  polar  array  that  varied  from  90  to  145  degrees, 
relative  to  the  upstream  jet  axis.  The  far-field  noise  was  measured  with  and  without  the 
microjet  injection.  As  can  be  seen  in  Figure  39,  there  was  little  or  no  effect  caused  by  the 
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microjet  injection.  In  fact,  there  was  even  a  slight  increase  in  the  overall  sound  pressure 
level  (OASPL)  at  both  140  and  145  degrees. 


Figure  38.  Setup  for  the  high  frequency 
forcing  experimental  setup. 


0  (degrees) 


Figure  39.  Plot  of  the  far-field  noise  radiation  for  the  high  frequency  forcing 
experiment. 


Even  though  the  results  suggest  that  high  frequency  forcing  did  not  produce  any 
favorable  results,  we  cannot  state  with  certainty  that  this  is  the  case.  There  are  a  number 
of  parameters  that  can  still  be  investigated,  such  as  the  size  and  number  of  microjets  in 
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the  nozzle.  It  would  also  be  beneficial  to  supply  each  microjet  with  its  own  high-speed 
valve  that  is  close  to  the  nozzle.  This  would  allow  for  individual  control  of  each  of  the 
microjets,  as  well  as  reducing  the  losses  that  undoubtedly  occurred  in  the  setup  used  for 
this  experiment.  Work  is  currently  ongoing  at  FMRL  to  develop  a  smaller  high-speed 
valve  that  can  be  used  for  these  types  of  applications. 


5.  SUMMARY 


A  novel  high-speed  jet  noise  suppression  technique  using  high-pressure  gas 
microjet  injection  at  the  nozzle  exit  was  developed  with  promising  results  using  the 
laboratory  scale  jet.  The  main  jet  parameters,  the  nozzle  pressure  ratio  and  the 
temperature  ratio  are  chosen  to  correspond  with  realistic  engine  operating  conditions. 
Keeping  in  mind  the  applicability  of  the  technique  to  full-scale  engines,  the  microjet  mass 
flow  was  kept  at  less  than  2%  of  the  primary  jet  mass  flow. 

One  of  the  dominant  sources  of  noise  from  a  high-speed  jet  is  intimately  related  to 
the  large-scale  vortical  structures  that  convect  at  supersonic  speeds  relative  to  the  ambient 
medium.  The  dominant  noise  from  this  source  consists  of  two  parts:  the  Mach  wave 
radiation  and  a  contribution  due  to  the  unsteady  flow  associated  with  eddying  motions 
whose  scale  is  commensurate  with  the  local  shear  layer  thickness.  The  noise  radiation  is 
most  intense  in  the  direction  given  by  a  uniquely  defined  convective  Mach  number  due  to 
Oertel5.  The  Mach  wave  radiation  that  is  quite  distinguishable  is  characterized  by  strong 
narrow  pressure  transients  in  the  far-field  microphone  signal,  referred  to  as  ‘crackle’,  first 
brought  to  notice  by  Ffowcs  Williams8.  It  is  quantified  by  the  skewness  factor  of  the 
pressure  time  signal  that  does  not  depend  on  the  scale  of  the  jet.  The  microjet  injection 
scheme  used  here  significantly  reduces  its  level  resulting  in  about  1.5  dB  reduction  in  the 
far  field  OASPL.  Additionally,  the  microjets  interfere  with  the  unsteady  vortical  flow  due 
to  large  eddies  so  as  to  reduce  the  mixing  noise  by  about  3  dB.  The  A-weighted  spectra 
that  incorporate  the  full-scale  nozzle  diameter  show  total  OASPL  reductions  of  about  6 
dBA  in  the  peak  radiation  angle.  The  microjets  are  effective  probably  because  they 
inhibit  the  formation  of  the  large  eddies,  the  consequence  of  basic  jet  instabilities,  which 
are  responsible  for  the  crackle  and  the  low  frequency  component  of  the  mixing  noise. 
When  the  jet  is  imperfectly  expanded,  the  commonly  observed  screech  tones  are 
completely  suppressed  by  the  microjet  injection. 

Another  major  component  of  the  jet  noise  is  due  to  chaotic  small  and  intermediate 
scale  turbulence  in  the  jet.  Although  it  is  difficult  to  suppress  this  noise  component,  our 
attempts  with  water  injection  using  microjets  appear  to  show  promise. 
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An  attempt  is  made  during  this  investigation  to  develop  Active  Noise  Control 
methods  for  cancellation  of  broadband  jet  noise.  A  demonstration  of  successful  control 
architecture  has  been  accomplished  using  the  colored  random  noise  cancellation.  During 
this  investigation,  an  attempt  was  also  made  to  develop  a  high  frequency  actuator  with 
based  Piezo-operated  valve  to  provide  flow  modulation.  The  valve  although  showed 
some  promise  failed  produce  sufficient  authority  to  effect  the  jet  flow  upstream  of  the 
nozzle  exit. 

Further  investigations  were  carried  out  to  design  and  built  an  actuator  using  the 
Hartman  resonance  tube  principle.  The  resulting  output  of  158  dB  from  these  actuators 
was  again  found  to  be  insufficient  to  modulate  the  flow  upstream  of  the  nozzle  exit.  A 
third  actuator  based  on  pulsed  microjet  appears  to  be  most  promising  for  the  current 
application.  A  specially  designed  micro  valve  operating  at  high  pressures  and  high 
frequencies  was  developed  and  tested  successfully.  The  valve  was  able  to  provide  a 
pulsed  jet  operating  at  nozzle  exit  mean  pressure  of  about  45  psia  at  2000  Hz.  Preliminary 
experiments  using  four  pulsed  jets  injected  upstream  of  the  nozzle  exit  in  an  open  loop 
control  configuration  showed  very  little  effect  on  the  far-field  noise.  However,  further 
experiments  are  being  conducted  under  a  continuation  ONR  grant  to  further  explore  the 
broadband  active  noise  control  of  jet  noise. 
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